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Introduction 


“I’m not happy with all the analyses that go with just the classical theory, because 
nature isn’t classical, dammit” 


This brusque statement was made by Richard Feynman in his keynote speech ‘Sim- 
ulating Physics with Computers’ in 1981 and heralded the new field of Quantum 
Computation. He stressed the fact that computers whose behaviour depend on quan- 
tum mechanical phenomena are perhaps more powerful than ordinary, ‘classical’ com- 
puters. This has proven to be true and has raised a growing interest among both 
physicists and computer scientists in the theory of Quantum Computing. 

In 1994 Peter Shor showed the existence of an algorithm for a quantum computer 
that can factor any number within polynomial time. This is a remarkable result 
because it is generally believed that on a classical computer this problem takes up 
an exponential amount of time, which is the reason that the difficulty of factoring 
is used in most modern cryptography protocols. More recently, Lov Grover demon- 
strated that a quantum computer can find an entry in an unsorted list of size N with 
a time-complexity proportional to only VN. It is not known how to do this for any 
computer as-we-know-it-today. 


We will apply the paradigm-shift from classical to quantum computation to the model 
of Cellular Automata. Cellular Automata are commonly used to describe discrete 
systems with a parallel and uniform time-evolution. ‘Conway’s Game of Life’ is a 
well-know example of such an automaton on a two dimensional grid. Just as Turing- 
machines are appropriate when considering sequential computation, cellular automata 
provide us with a theoretical abstraction of massive parallel computation. Cellular 
automata are also used in physics, biology and other areas to describe systems such 
as fluids, gases and DNA-sequences. 

After investigating some of the typical characteristics of this model, it will be 
shown that there exists a Universal Quantum Cellular Automaton that can simulate 
any other automaton with only linear slowdown. 'This result may be of use for the 
actual construction of a quantum computer because several authors have suggested 
that it may be easier to construct a Quantum Cellular Automaton than a Quantum 
Turing-machine. 


Chapter 1 


Quantum Physics 


We shortly describe some concepts and laws of quantum mechanics. We will re- 
strict ourselves to the fundamentals which are necessary and sufficient to understand 
quantum computation. 


1.1 The Two-Slit Experiment 


As a start we will look at Young's two-slit experiment which shows us the basic properties of quan- 
tum physics (Figure 1.1). After describing this well-known experiment in words, the bracket no- 
tation will be introduced. The mathematical background of this notational tool will be explained. 
We refer to the standard introductory books for a more thorough and detailed explanation of 
quantum physics [10, 14, 24]. 

The set-up of the two-slit experiment is as follows. A source S sends a photon through a wall 
with two-slits A and B on a projection screen P. The probability of a photon arriving at position 
x on this screen will depend on the value of x. When we repeat this experiment for a large number 
of photons this probability-distribution on P will produce an interference pattern which can be 
explained by the wave-like properties of light. 

The problem is that this interference also occurs when the source emits photons at such a low 
rate that only one photon at the time travels from S to P. This means that a photon can interfere 
with itself. To understand this phenomenon we have to look more precise to what is happening. 

For every photon that arrives at x there are two possibilities for its history: it either went 
through slit A or slit B (from the source S). After the photon has been observed at x it is 
impossible to determine if it used A or B to pass the wall. We therefore say that the two routes 
“from S, through A, to x” and “from S, through B, to x" are indistinguishable. This is a necessary 
feature of the experiment because when we do track down the way the photon goes, the interference 
pattern vanishes. (This can be done by setting up two detectors at the two slits.) In short: a 
photon that goes through specifically one slit cannot interfere with itself. 


1.2 The Bra-ket notation 


We now know that the interference depends on the set of possible paths a photon can travel from 
S to x. This can be described by the bra-ket notation. This is the conventional way of describing 
quantum mechanical events and was introduced by Paul Dirac [17]. The mathematical basis of 
this formalism will be explained in the next section. 

When it is possible to go from a configuration X to configuration Y we can describe this with 
a non-zero probability amplitude (Y |X). This amplitude is a complex number with a norm < 1. 
The (Y | part of this expression is called a bra, the | X) part is a ket. The whole term is therefore 
called a bra-ket. Inside this bracket are descriptions (note: not mathematical expressions). A 
bracket should be read from right to left to understand its meaning. If it is not possible to go 
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projection screen 
wall 

x 

source A 
S 
B 
interference pattern 
P 


Figure 1.1: Young's two-slit experiment. The source S transmits photons through a wall with 
two slits (A and B) on a projection screen P. The arrival of a single photon at the screen 
P is described as the superposition of the two possible trajectories allowed by the two slits. 
Because the probability amplitude of arriving at position x via A differs (in general) from the 
amplitude of arriving via B, an interference pattern will appear. If a measurement is performed 
in order to determine through which slit the photon has traveled, the interference pattern at P 
will disappear. This is because the two trajectories are now made distinguishable which destroys 
the initial superposition of the photon. In reality distance between the slits A and B has to be 
much smaller to get a noticeable interference pattern. 


from X to Y the amplitude (Y|X) will be 0. We say that X and Y are mutual orthogonal when 
(Y|X) =0. 

By multiplying the amplitude (Z|Y) with the amplitude (Y|X), we get the amplitude of going 
from situation X to Z, thus: (Z|X) = (Z|Y) (Y|X). 

A probability amplitude tells us the probability | (Y |X} |? of measuring Y given the initial 
configuration X. The amplitude (Y| X) cannot be measured directly. 

Because the labels X and Y describe situations they can sometimes be decomposed into 
more specific descriptions which are mutual exclusive. If we take the example for which |X) — 
|X, or X5) then the amplitude of going from state | X) to state (Y| can also be decomposed. If 
the parts | X1) and |X2) are indistinguishable we can add the amplitudes such that: 


YIX) = (YX o Xə = (YMIXO-c-IX39) = (YIX) + (Y1X2) 
When the parts are distinguishable we cannot add the amplitudes. The same reasoning goes for 


the bra, for example: (Y |X} = (Y1 or Yo|X) = (Y1| X) + (Yo|X). 
Let us now restate the two-slit experiment with the use of this notation. 


EXAMPLE 1.1 The amplitude of an emitted photon to reach position x at the screen is described 
by: 


( arrives at x | emitted from S) =  (x|S) 


As we have seen this trajectory can be divided into two indistinguishable routes “going through 
slit A” and “going through B”, therefore: 


(x|S') = («| from A) ( through A |S) + (x| from B ) ( through B |S} 
= (a|A) (ALS) + (B) (BIS) 
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The probability of a photon to reach x from S is thus calculated by: 
| (|S)? = | (aA) (AIS) + (2| B) (BIS) |? 


This shows how the amplitudes (x| A) (A|S) and (x|B) (B|S) can interact and thereby increase or 
decrease the value | (x|S) |. If the A and B-part have the same phase, the probability is amplified. 
Take for example the values: (x|A) (A[S) = (x| B) (B|S) = 1/2. The probability of a photon to 
reach x via A or B individually equals |1/2|? = 1/4. But the total probability is |1/2 + 1/2|? = 1. 

The two parts ‘cancel’ each other out if they have opposite phases. To see this we have to 
change (r|B) (B|S) into —1/2 such that | (r|S) ^ = 0. This phase-dependent behavior explains 
the interference pattern in our experiment. The phase of a photon to go from a slit to a position 
x on the screen depends on the distance between the slit and x. Because the distances A — x and 
B — x are different for every x, the two amplitudes (z|A) and (x|B) sum up differently for every 
position 2x. 

When we make the A and B parts distinguishable (for example by using a photon-detector 
at the slits) this expression changes. Because we cannot add the amplitudes (r|A) (A|S) and 
(|B) (B|S) anymore, the probability of measuring a photon at x now equals: 


l(x|S) ^ = |( arrives at x through A |S) |? + | ( arrives at z through B |S) |? 
| (214) (ALS) |? + | (zB) (BIS) |? 


We simply have to add the two probabilities. 
The difference between the two calculations lies in the fact that amplitudes are complex-valued 
numbers but probabilities are always in the domain [0,1] C R. © 


1.3 The Mathematics of Quantum Physics 


We will now make the bracket notation more precise and meaningful in a mathematical sense. 
This can be done by defining a complex valued vector space such that the |-) and |-) expressions 
are vectors with an inner product (-|-). For our purposes it will be sufficient to use the following 
two definitions of a Hilbertian space and a state space. 


DEFINITION 1.1 (Hilbertian space) For every set of basis states D, the Hilbertian space (9(B) 
is the complex-valued linear space on the domain B with a bounded norm. In other words, for 
every 


X = Mot 
EEB 


with og € C, the vector X is an element of /5(B) if and only if (a* is the complex conjugate of 
a € C) 


X agag < oo 


EEB 


This space is equipped with an inner product (-,-) : (9(B) x £9(B) — C and a norm ||-|| : €9(B) —> R 
defined by: 


(X,Y) = (Zas ; Yao) = Mie 
EEB EEB EEB 


and ||X|| = y (X, X) for every X,Y € ¢2(B). By definition all the vectors in B are mutually 
orthogonal and have norm 1. 
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Note that we do not claim that this is the formal definition of a Hilbert space. The above described 
£5(-) space is only an example of a Hilbert-space. It can be shown that the inner product and 
norm have the properties: 


(X, aY + 8Z) =a(X,Y) +8 (X,Y) |X|] = 0 

(X,Y) = (Y, X)" || X || = 0 if and only if X = 0 
(X, X) 20 IXI + IYI 2 X + YI 

(X, X) = 0 if and only if X 20 llo-X || = Ja] || XI| 


for every o, 8 € C and X,Y, Z € (9(B). Which are exactly the properties we want them to have. 
'The complex-values in combination with the notation of the inner product already forecasts 
the next definition which formalizes the bra-ket notation. 


DEFINITION 1.2 (State space) Given a set of basis states D, the state space Hg equals the 
Hilbertian space ¢2(B) restricted to the vectors with norm 1. The bracket notation enables us to 
describe vectors by (X| and |Y} € Hp. If we restrict the inner product on (B) to the domain 
Hpg, we define with this notation: 


(xY) = (X|) 


The Hilbertian space £5(B) is a complex, linear space spanned by the basis set B. Every vector 
is in £2(B) is therefore a linear combination of basis states € B. Because the state space Hpg is 
a subset of /3(B), every state |X) € Hp is also described by a linear combination of basis states. 
A state does not determine its basis set. It is therefore possible to have two different basis sets 
A £z B with Hg = Hg. 


1.4 The Superposition of Basis States 


The possible configurations of a physical system are described by the state space Hg which is 
spanned by a set of basis states B. The number of basis states (which equals the dimension of 
Hp) can be infinite and even uncountable (like the arriving of a photon on a position z € R in our 
example). From now on we will only look at finite dimensional state spaces. A basis state £ € B 
will be denoted by |£) to preserve the bracket notation. 

As a consequence of this, every state |.X) can be described by a linear combination on the basis 
states B: 


IX) = led 


EB 


with a € C. It is said that X is in a superposition of basis states B. 
Because B is an orthogonal basis set with all vectors norm 1 we have for every £, x € B: 


O 1 ify =€ 
&b = li size 
By multiplying the linear combination of X on the left with a basis state (y| we therefore get: 


(IX) = je) = ay 


EEB 


for every x € B. This shows us how the inner product in the Hilbertian space 2(B) relates to the 
probability amplitudes in quantum mechanics and vice versa. By definition is holds that 


IX) = Sole) (EX) 


EEB 
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for every state |X) in the state space Hg. The probability of measuring a state |X) in the basis 
state £ equals | (£| X) |. By the ‘norm 1’ constraint on the state space Hg we see that the overall 
probability of measuring a configuration X in a basis state is 1: 


$e = 1 


EEB 


for every |X) € Hpg. This is the normalization restriction which will play an eminent role in this 
thesis. Only states and transformations which respect this restriction are sensible from a physical 
point of view. We say that a state space and its transformations have to be proper or well-formed. 
This normalization condition for Hp is a restriction on the Hilbertian space /2(B). The classical 
state space B is a subset of Hg where every amplitude has the value 0 or 1. We therefore can 
write: 


B Ç Hp C Q(B) 


1.5 Evolution of Quantum Mechanical Systems 


If we want to be more precise about the (time)-evolution a system undergoes, we have to use 
transformations which describe this evolution in the state space. The transformations are con- 
ventionally called operators or time operators. If under the influence of a operator U a systems 
evolves from state | X) to a state |Y), we denote this by: 


UX) = JY) or (Y|UX) = 1 or X” —v quj 


If we want this operator to obey the laws of quantum mechanics U has to be both linear and norm 
preserving. It has to respect the normalization condition and the superposition principle of the 
state space Hg. We therefore can write 


UX) = Ul Scale) | = ie Ulo 


EEB EEB 


for every |X) € Hg. Because U is a linear transformation it can be described by a matrix My. 
If we want this matrix to be norm preserving it has to be a unitary matrix. In order to define 
unitarity we first have to define the hermitian conjugate of a matrix. 


DEFINITION 1.3 (Hermitian conjugate) The hermitian conjugate of a matrix M is the matrix 
M! with [M`]; = [M]}; for every index i and j. 


This enables us to define: 


DEFINITION 1.4 (Unitary matrix) A unitary matrix M is a complex valued matrix such that 
the hermitian conjugate M' is the inverse of M and vice versa. Therefore: 


M.M'! = Mi-M = I 


A unitary transformation is thus defined by a unitary matrix. From now on we will make no 
distinction between a function U and the corresponding matrix U. 


DEFINITION 1.5 (Unitary transformation) A transformation U : Hg — Hpg of the state 
space Hg is unitary if and only if it corresponds to a unitary matrix U € C?*4 with d the 
dimension of Hp. 


Because every unitary matrix U has an inverse U' it defines a bijective or reversible transformation. 
The classical unitary transformations are a subset of the unitary transformations and are described 
by the unitary matrices U with [U];; € {0,1} for every i,j. See the appendix for the definition of 
a unitary transformation on a infinite dimensional state space. 
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1.6 Expanding the State Space 


Consider two state spaces Hy and Hw spanned by V = (v1,... , Un} and W = {w1,... , Wm}. If 
we join these spaces we get an expanded space which will is expressed as the tensor product ® of 
Hy and Hw. This gives us a new state space Hz with dimension nm: 


Hz = H(vxw) = Hv @Hw 


This space is spanned by the basis vectors zi; = v; ® wj. In the bracket notation this is usually 
denoted as: 


|v) lw = lw) = lea) 


for every v; € V, wj € W and zi; € Z. The tensor product & is a linear function which determines 
all the vectors of Hz by: 


Eo |X)= Cali) and  |Y)-Y]8;w;) 
i=l j=1 


then  |X)&|Y) = JM aiGjlvi,w;) 
i,j=1 
In vector notation this can be visualized as: 
ay fy 
oi|Y) a B» 
a2|Y) s 
|X)e|Y) = : = andi EC 
On|Y) l 
On Bm 


If the states | X) and |Y} are known states, we will denote the joint state | X) & |Y ) also by | X, Y). 
This is called a product state. (Because X and Y are just labels this is not a mathematical rule, 
just a notational custom.) The set of product states of a state space Hy & Hw = Hz is a proper 
subset of Hz. The states which are not product states are called entangled states because there 
is correlation between the V-part and the W-part of the vector. 

If a state space Hz can be decomposed into two (or more) spaces Hy and Hw such that 
Hz = Hy & Hw, we refer to Hz as a product space. Now V and W define the ‘subspaces’ or 
‘subsystems’ of Hz. We will encounter these subsystems later on, when we are discussing the 
notion of entanglement. 

If P is an operator on Hy and Q is an operator on Hw, then their combined behavior in Hz 
is also a tensor product: 


P|X)eQ|Y) = (PeqQ)X.vY) 
Because P and Q are matrices, the tensor product P & Q can be calculated by: 


PQ PQ > Png 
puQ P2Q +--+ PmQ 


€ (uomo 


PEQ = 
PriQ Pn2Q > Pm 
The most important properties of the tensor product are (assuming that the dimensions of P,Q, R 
and S are compatible): 
If u € C then (uP) & Q = P & (uQ) = u(P & Q) P®(QSR)=(PeQ)@R 
(P+Q)@®R=P@R+QOR (P ®Q)(R®S) = (PR QS) 
PS(QrR)-PSQM-PGR (PEQ) = (Pİ 9 QÏ) 
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Tensor products for matrices are also called “direct products” or “Kronecker products” [25, 33, 47]. 
The tensor product should not be confused with the Cartesian product “x”. A tensor product of 
two spaces H4 and Hp is induced by the Cartesian product of the two basis sets A and B. 


1.6.4 Notational Customs 


The following notational customs will be used in this paper. When necessary, a brief explanation 
will be given. 


repeated tensor products By P"! (with n € N+) we mean P & P®---® P (n times) and by 
definition P] = 1 € C. Because a tensor product does not commute, we have to define 
explicitly: 


i=r—1 
Q Pi = Q P = M*8ƏP 8:8 Pı 
i€Z, i=0 


for r € Nt, and because P = 1 also: 


QP: = 1 eC 


state space If we use the basis set B, the corresponding state space will be denoted by Hpg. 


basis set/states The set B is will be treated as a subset of Hg. The basis states in B will also 
be named canonical states. 


ket notation For every basis state £ € B we will write |£) when £ is used in the context of state 
space Hpg. 


long ket notation Because |x) & |y) implies that the x and y part are not entangled, it is 
erroneous to state |r, y) = |x) & |y). In general -if entanglement is involved- we have to 
describe the state by one large ket. To avoid cumbersome expressions, the following notation 
is introduced: 


leas 4%) = Je.) = de) 


1.7 The Effects of Measurement 


A quantum mechanical system is fully described by its state vector. This state however, cannot 
be measured directly. That is, we cannot observe the amplitudes of the different base vectors. 
Only some specific information of this state can be observed. These observables are characterized 
in the following way. 

Consider a state space Hz with a state |W) € Hz. An observable corresponds to a set of 
subspaces 7tw,,w,,..., Hw, € Hz with: 


e for every i, j, if i Aj then Hw, 1 Hw, 
e for every state |V) c Hz: |W) = a o;v; with |w;) € Hy, for every i 


An observation along these subspaces will give us a result ‘j’ which corresponds to a single subspace 
Hw. This obeys the following rules: 


e The probability of measuring j equals [o;;|?. 


e If j is measured, the state collapses according to |Y) ~~ |v;) (the new state vector is 
normalized again) 
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Because of the second rule all the information about the amplitudes a is lost. This shows that mea- 
suring a quantum-mechanical system is in general an irreversible process because the superposition 
is disturbed (hence the “~»” symbol). 


EXAMPLE 1.2 Take a two spin system Hz with Z = {|ff),|TL), L1), |LD)) with a state |W) in 
the superposition: |Y) = $|11) + 2111) + $111). A measurement on the first spin divides Hz into 
the two subspaces Hw, and Hy,, with Wi = (11),|11)) and Wo = {|IT), |ll)}. The Wi space 
corresponds to measuring “up” and Wọ with measuring “down”. 

Along these subspaces, the state vector is described by: 


W) = {Fin malty b+ Stuy} = V5 a) + le) 


with |W) € Hw, and |v2) € Hw,. As a result, the outcome of the measurement has two possibil- 
ities: 


1. first spin is “up” with probability 5; the new state vector becomes: |W’) = 4V5{|11) +2|T|)} 
2. first spin is “down” with probability 4; the new state vector becomes: |W’) = |||) 
© 


This example shows that if a state |X) € HA8&Hp has entanglement between the two subsystems, 
an observation on the A-part will also influence the B-part of the state. 

Our first lemma tells us how to calculate the amplitudes of a subsystem Hg in an expanded 
state space Han. 


LEMMA 1.1 Given an expanded state space Hon whose states are described by |x1,... , v4) € 
Han. For every X € Han, the probability amplitude to have |Y) € Hg on the xı position of the 
state | X) equals: 


(a -Y|X) = SO im =Y,22n = Z|X) 
ZeQn-! 


PROOF. Let Q’ be a basis set with Hg = Hgo and Y € Q'. Because Q’ is an orthogonal set, we 
have for every q € Q*: 


7 7 | f1 ifq=Y 


The state space Hg. gn-1 equals the state space Ho» and therefore (by summation over the new 
basis set Q’ x Q"-1): 


(n-Y[X) = $7 X5 (m1 =Y]|tı =q, tzn = Z) (£1 = q, tom = Z|X) 
qEQ' zeQn-1 
= 5 (£1 = Y, Gon = Z| X) 
ZEQr-t 


Chapter 2 


Quantum Computing 


In this chapter we will describe the basics of quantum computation. This is done by 
defining quantum bits, quantum registers and quantum gates. The model of quantum 
Turing-machines will also be mentioned. 


2.1 Quantum Memory 


2.1.1 Quantum Bits and Quantum Registers 


Consider a two-state system B for which we have labeled the basis states “0” and “1”. The state 
set Hp is therefore be described by the configurations 


|X) = aol0) + o4 |1) 


for every ao, a, € C with |ao|?+]a1|? = 1. Such a state |X) will be called a quantum-bit or qubit. 
Because the canonical states |0) and |1) will have their natural meaning, this systems is usually 
indicated by H{0,1}- 

If we expand this state space to that of a system whose basis set is described by {0,1}” we 
get the definition of a n-qubit system. The possible configurations of such a quantum register are 
covered by the expressions 


IX) = M, od? 
&cto.1)» 


which obey the normalization restriction. The state space of an n-qubit system equals the tensor 
product of n separate qubit systems: 


[n] 


Hoar = Hon 871 2 @Aoy = Hion 
$$ 
n times 


Qubits and quantum registers are used to describe the memory of quantum computers. The canon- 
ical configurations of an n-qubit system are the ‘classical configurations’ {0,1}". The canonical 
configurations are a basis set for the system 7Hy9,1}». 


2.1.2 Entanglement 


When the amplitudes of an n-qubit configuration define a correlation between the individual qubits, 
we say that the qubits are entangled. This correlation can effect the outcome of measurements 
performed on the system in a way that can not be understood in a classical sense. 
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EXAMPLE 2.1 Take the following configuration of a two-qubit system: 
1 
|) = 5v24]0,0) + |1,1)} 


When we measure the value of the first qubit, the wave function will collapse to one of the two 
possible states: 


| Wo) = |0, 0) or |V) = |1,1) 


This means that the value of the second bit is determined by the outcome of our measurement 
on the first bit: there is a correlation between the bits. This is not always the case. If the initial 
two-qubit system is a tensor product of two individual qubits, then a measurement on the first bit 
does not affect the second bit. For example: 


1 1 
|) = 5{l0,0) + 10,1) + 12,0) +1,1} = 51000-1101 8 (02 + |12)} 
An observation on the first bit changes the system into one of the configurations: 
1 1 
|Yo) = |01) & 5V 24102) +|12)} or |Wi) = |1) ® 5 V2 02) +|12)} 


which does not affect the second bit. © 


The notion of entanglement is important because it illustrates the influence of a measurement 
on the system and is typical for quantum-mechanical systems. It also enables the occurrence of 
interference in a quantum computer which defines the fundamental difference between probabilistic 
and quantum computing. 

Entanglement can not always be seen by the probability of measuring certain basis states. As 
a counterexample: the two qubits in the following state are entangled: 


I) = 5 {10,0) + 0, 1) +11,0) - [1,1)} 


although there is no ‘correlation’ to be found when measuring the bits. 


2.2 Quantum Gates 


We will now look at gates which operate on qubits: quantum gates. In the previous chapter we 
discussed unitary transformations operating on a state space. A quantum gate is a system that 
performs such a proper transformation on a register of qubits. Every proper quantum gate that 
operates on a n-qubit system is therefore described by a unitary matrix M € C4*4 (with d = 2” 
the dimension of the state space H1o,1}»). An example of a gate that operates on one qubit is the 
NOT-gate. 


EXAMPLE 2.2 The NOT-gate is described by the unitary matrix which operates on the two 
dimensional state space Hgo} with 


(Dem m (2) 


The matrix Myor therefore equals 


|1) 
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The ‘traditional’ behavior NOT(0) = 1 and NOT(1) = 0 is shown by the matrix multiplications: 


Ment) = (9 5)-(9) = (1) = m 
Meb = (2 5)-(9) = (5) = v 


In general we have for this gate: 


Motos) ul - (1 a (82 )=( 88) = eo + aolt 


and 


ay 


with o9, o1 € C and laol? + lail? — 1. o 


'The last equation in this example shows that a quantum gate has to respect the possible super- 
position of a quantum register. This enables us to define gates whose behavior is impossible with 
classical gates. An example of such a true quantum gate is the following: 


EXAMPLE 2.3 Define a one-qubit gate Q by: 
lfí1-i 1+i 
LION ( 1+i 1-i ) 
When applying this gate to the canonical states |0) and |1), we get: (A= $ — 4; |A|? = 3) 


Q|0) = lO) + A*|) and QI1) = A*|0) +.A}1) 


In both cases the result is a perfect ‘fifty-fifty’ mixture of |0} and |1) (there is only a phase difference 
between Q|0) and Q|1)). If we apply the Q-gate a second time on this result we get: 


Q{Q|0)} —|1) and Q{Q|1)} = |0) 


which is the NoT-function again. This is affirmed by the matrix of Q?: 


2 0 1 
e - (io 
We therefore say that Q is the “square root of NOT”. Although the input and output of Q? 


‘classical’ it is not possible to have a classical gate Q with the same behavior. This is proven by: 
det(Q)? = det(NoT) = —1 and thus det(Q) € (i, —i) which is impossible for a classical gate. © 


'The reversible gates with classical behavior are a proper subset of all possible quantum gates. 
Non-reversible gates are not proper because the corresponding matrix has to be an invertible 
matrix (M - M! — 1). 

2.2.1 Proper Quantum Gates 

It is not always obvious if a gate is well-formed or not. We therefore will use the following ‘rule of 
thumb' to certify that a gate with a certain requested behavior is proper. 

DEFINITION 2.1 (Requested behavior) The requested behavior of a gate M is described by a 
finite list of input and output values ((x1,y1).... , (£k, yk)) such that M (x;) = yi for every i. 


'This definition is used in the next lemma about proper quantum gates. 
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LEMMA 2.1 If the requested behavior of a gate M : Hga — Hp» is described by a list 


(21,91), CE (£k, Yk)) 


with x; L x; and yi L yj for every i # j, then there exists a proper quantum gate with M(x;) = yi 
for every 1 € i € k. 

PROOF. Let the set of basis states be denoted by B" = {&,...}. Because z1,... , zy are mutually 
orthonormal, there exists a unitary transformation P : Hga — Hen with P(£;) = x;, for every 
1<i<k. For the same reason there also exists a unitary R such that R(£;) = y;. If we define 
M = R- P', M will also be unitary, with for every 1 < i < k: M(z;) = R(P'(z;)) = R(&) = yi. 


2.3 Quantum Gate Circuits 


An acyclic circuit of quantum gates is called a quantum (gate) circuit or quantum gate array. 
'The general behavior of a quantum circuit can be calculated with the matrix formalism in a 
straightforward way. 

Consider two gates A and B that operate on two disjoint subsystems of n and m bits: 


A: Hoy = Hyo,1yn and B: Hom — Hom 


The combination of A and B defines a transformation on the state space 719,1» € 79,» = 
Hyo,1}n+m described by the tensor product A & B 


(A& B): oim — Hoi 
with 
(A@ B)X,Y) — AX) 6 B|Y) 


for every X € (0, 11* and Y € {0,1}™. 

When two layers of gates act as a sequence on an n-qubit system, the general transformation is 
calculated by the product of the two defining matrices. This rule was already used in the example 
of the square root of NOT gate. To summarize: 


1. When two gates are sequential, multiply the matrices 
2. When two gates are parallel, take the tensor product of the matrices. 


See Figure 2.1 for an example of a small quantum circuit. The characteristics of the quantum 
circuit model were first investigated in articles by David Deutsch [16] and Andrew Chi-Chih Yao 
[60]. 


2.3.1 A Universal Quantum Gate 


Because quantum arrays define unitary transformations on qubit systems, it is natural to ask 
^What kind of quantum gates do we need to implement an arbitrary unitary transformation?" 
For classical non-reversible computation we know that it is sufficient to have AND, OR and NOT- 
gates. The quantum case is a more difficult problem to solve. 

It is impossible to construct every unitary transformation exactly with a finite set of basic gates 
because the number of possible transformations is uncountable. We therefore have to approximate 
the transformation within a bounded error. This means that given a transformation U and an 
arbitrary small number £ > 0 we can make a circuit U’ which simulates U within the allowed 
error £, that is |U — U'|| < &. If there exists a quantum gate by which we can construct any 
transformation within a bounded error, we will call this gate a universal gate. 
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layer with gates: A,B,C 


layer with gates: D,E 


Figure 2.1: A small quantum gate circuit with two layers and five gates. The action of the first 
layer is described by the tensor product A & B & C. The second layer equals: D & E. The 
transformation defined by this circuit is therefore described by: (D & E)-(A& B&C). 


It was proven by David Deutsch that there exists a universal gate that operates on three bits 
[16]. After that, Adriano Barenco[3] and David DiVincenzo[18] showed the existence of a two bit 
gate which is universal in is computational power. For a detailed expose of this subject the reader 
is best referred to the famous ‘article by nine authors’ by Adriano Barenco et al. The quantum 
circuit model has proven to be the most convenient model to study quantum computing. 


2.4 Quantum Turing Machines 


'The quantum Turing machine provides us with an alternative model for quantum computing. 
'There is strong relationship between the characteristics of such QTMs and reversible Turing ma- 
chines [5, 15]. The following definition has its origin by Ethan Bernstein and Umesh Vazirani 
[6]. 


DEFINITION 2.2 (Quantum Turing machine) A Quantum Turing Machine M is defined by 
the tuple (K,3,0), with X a finite alphabet set (including a blank symbol O), K a finite internal 
state set and ô : K x x (—,] x K x X 2 C the finite state control. 


The set of basis states is described by the states |Y) = |k, h, V) € L(K x Nx Q*) 2 S, where k 
indicates the internal state of the QTM; h defines the head-position on the one-way infinite tape 
and « describes the tape-content (with only a finite number of non-blank symbols). The global 
transition function M : S — S now obeys (for every k, h, v), q, h', k' and y’): 


(K', h— 1, Vo:h—1 x qx Wn+1:...\M|k, h, v) ó(k, Wh, =,q, k^) 


(K', h +1, o1 xq Xx Wn41:...\M|k, h, v) = O(k, Yh, —,q, k^) 
(Ku! |M|k, hv) = 0 otherwise 


(Note that by each computational step the head must move.) By superposition on the set of 
basis states, M describes a transformation on the Hilbertian space S. If M describes a unitary 
transformation, the QTM M is said to be well formed or proper. 


NOTATION 2.1 We will use the following notational conventions: (|X),|Y) € S and M a well 
formed QTM): 


e Instead of |k, h, v) we may also write |k, o, ... , Vni, Yh, has ++): 
e By |X) — y |Y) we mean: M|X) = |Y} 


e If there exists a t € N with M'|X) = |Y) we denote this by |X) =m |Y) 
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e The internal state is often described by a function-value ' ...'. For example, with ' start ' 
and ' halt! we mean the starting and halting state € K of the QTM. This is done without 
further specification. 


2.4.1 Well Formed QTMs 


Because QTMs will not play an important role in this paper, we will only summarize the main 
characteristics of quantum Turing machines. First of all we have to take a closer look at the well 
formedness issue. To handle this restriction on the infinite dimensional Hilbertian space S, Ethan 
Bernstein and Umesh Vazirani translated it into the following three constraints. 


LEMMA 2.2 A quantum Turing machine is well formed if and only if the finite state control ó 
obeys the following local conditions (where » ^, p |K, h, v») describes the summation on the basis 
states of S): 


1. For every basis state X, the transformation M is norm preserving: 


XO kh yM = 1 


kh 


2. If X and Y are different basis states with identical head positions, then M(X) and M(Y) 
are orthogonal: 


So IMIE h, 4%)(k,h,Y|M|X}) = 0 
k,h,p 


3. For every basis state X with head position h the *—' and ‘>’-part are mutual orthogonal: 


X (k,h+1,Y|M|X)(X|M|k,h- 1,4) = 0 
kw 


PROOF. See the original paper [6]. 


A QTM with a deterministic head position always satisfies the third constraint. We will use this 
in the following lemma which shows us that any finite dimensional unitary transformation can be 
simulated by a quantum Turing machine. 


LEMMA 2.3 For every unitary transformation M on ?H(o,1» there exists a well formed QTM T 
such that for every |) € Hyo,1}» we have: 


|"start7,0) @ |b) >r |Thalt7,0) & M|v) 


PROOF. We will restrict the definition to the values v € {0,1}” (by superposition this will impose 
the desired transformation on the whole state space H{o,1}» )- Because there is only a finite number 
(2”) of possibilities, we will use an exhaustive look-up table. First we will read the value from left 
to right (thereby simultaneously erasing the tape contents). After this part, the QTM ‘splits’ into 
the desired amplitudes and writes the various basis states on the tape. After a small ‘shuffle’ (the 
head is not allowed to stay stationary), the machine halts. With me, = (£|M |v) this is formally 
described by the sequence: 


\start ^, 0) & |v) — [vo 1,1) 8 |0, v1. 1) 
First we replace the n bits of the input string w to a corresponding state of the QTM ... 


> r|, n) & |0,... ,0) 
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after which the output value is calculated by a finite look-up table. This produces a superposition 
on the internal state of the QTM with amplitudes m. 


TT 5 Mey| v,£ n — 1) 8 |0,... ,0) 
&c(0,1)» 
By moving the head backwards, this output string is written on the tape of the Turing machine. 
MP yo Mey| Y, Eo:in-2 1,2 — 2) @ (Orisa yos) 
&c(0,1)» 
By doing so, the superposition of the machine T is transported to a superposition on the tape, 
>r 5 mey | W, £o ,0) & [0,61 1) 
&c(0,1)" 
Finally, the QTM is returns to one basis state ‘halting ... ’ which shuffles the head one last time. 
—r M; |Thalting ... 7,1) & meyl£) 


fe{0,1}" 
—r |' halt , 0) & M|v) 


With lemma 2.2 and the unitarity of M we can verify the well formedness of this QTM. Because 
M is norm preserving, the m values will also be proper (requirement 1). The ^... -function will 
be defined in such a way that v 4 wv' implies |47} L | w^. The transformation M is angle 
preserving, and therefore M|) L My’) which certifies the second restriction. Because this QTM 
has deterministic head position, the third requirement is trivial. 


It was shown by Ethan Bernstein and Umesh Vazirani[6] that there exists a Universal QTM that 
can simulate any other QTM with only polynomial slowdown. The equivalence of the quantum 
circuit model and QTM-model was proven by Andrew Chi-Chih Yao[60]. 


2.5 Algorithms for Quantum Computers 


What are quantum computers good for? When we want to simulate a quantum mechanical system 
a quantum computer could be very useful. On a classical computer we would have to compute 
the evolution for every basis state after which the final result is obtained by a summation of the 
calculated amplitudes for each basis state. For an n dimensional state space this calculation has a 
time complexity proportional to n. Because the dimension of the state space grows exponentially 
in the size of the system, this is a time consuming procedure. With a quantum computer this 
problem of large state spaces is solved by allowing the quantum computer to enter an equally sized 
superposition. By doing this we do not have to do the calculation for every basis state because 
the computer does this in parallel. This ‘quantum parallelism’ is also the main idea behind the 
algorithms that have been constructed to solve traditional computational problems in a faster way 
than is possible (or known to be possible) on a classical computer. 


EXAMPLE 2.4 Assume a function f : {0,1}" — {0,1} and a quantum algorithm M with 


M|£,0) = |&f(&) 


for every £ € {0,1}". If we provide this algorithm with an equally distributed superposition 
of input states, we can calculate the superposition of outputs with the same time and space 
complexity by: 


M| M z7^7k0| = Y) 27^|,f(6) 


EE{0,1}” €€{0,1}” 
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A measurement on this superposition gives us one of the basis states |€, f(€)) with a probability of 
2^". Because we cannot control the probabilities of a measurement it would be misleading to say 
that we really ‘know’ all the possible outcomes of a function. In general it is impossible to force 
the superposition into a basis state with f(£) = 1. In order to take advantage of this superposition 
we have to use the interference phenomenon on the amplitudes of the basis states [7, 22, 30, 51]. 
If this is possible and how to do this depends on the function f we have calculated. o 


In 1994 Peter Shor showed how a quantum computer could be used to factor numbers with only 
polynomial time complexity [49, 50]. Not only does this algorithm decide if a number is prime or 
not but it will also give the prime factors of a composite number. It is generally believed but not 
known that on a classical computer this problem requires an exponential amount of time to solve. 
Moreover, because of this intractability, factorization is commonly used in modern encryption 
schemes such as RSA. 

Another result was achieved more recently by Lov Grover [26, 9]. He proved that for any 
function f : Zn — {0,1} it possible to find a number i with f(z) = 1 (if such a number exists) 
with a time complexity proportional to yn. This quantum searching algorithm does not assume 
any knowledge about the function f. Although this is not an exponential speed-up it is likely that 
every deterministic or probabilistic algorithm will have a time complexity proportional to m for 
this task. 


Chapter 3 


Classical Cellular Automata 


In this section we will give a brief description of cellular automata and its character- 
istics. Special attention will be paid to the subclass of reversible cellular automata 
which play an important role in the the theory of quantum cellular automata. 


3.1 Cellular Automata 


The model of cellular automata (CA) is used to describe the behavior of discrete systems with 
a uniform and parallel space/time behavior. Given a space S and a local state set Q, a CA will 
describe a function Fs : Q9 — QÙ which defines the time-evolution of an initial configuration 
X € Q?. Consequently, the configuration at time t € N is described by the expression F'( X). 

Before giving a formal definition we will first look at a typical example of a one-dimensional 
cellular automaton to make ourselves familiar with some important characteristics. 


EXAMPLE 3.1 Consider a one-dimensional, two-state CA F such that the successor X' of a 
configuration X € (0,117 is determined by a local function f according to: 


X! = F(X) 


(gx) 


icZ 


II 


= Q) f(%i-1, Xi Xa) 


icZ 


This local function f : Q? — Q is defined by a finite table: 


f(0,0,0) = 0 f(0,0,1) = 1 
f(0,1,0) = 0 f(0,1,1) = 1 
f(1,0,1) = 0 f(1,0,0) = 1 
f(1,1,1 = 0 f(1,1,0) = 1 


which implements the addition-modulo-two-rule: f(£x—1, £o, £1) = x1 9 21. 

If we unfold the time parameter t of the global function F$ as an additional space-dimension, 
we get a (1 + 1)-dimensional structure which shows us the space/time behavior Z x N  Q of the 
CA on an initial configuration X: 

(i,t) e [Fix], 
Figure 3.1 shows us the result of this transformation. With the convention that time is going 
downwards each time-step can be identified as a layer in such a structure. 
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cell values at time t 


neighborhood scheme 


local functions 


cell values at time t+1 


neighborhood scheme 


local functions 


Figure 3.1: A part of the automaton in Example 3.1. The local states X; at time t are the 
input values for the local functions f which determine the next configuration at time t+ 1. The 
mapping between the cell-values and the local functions is described by the neighborhood scheme 
of the automaton. 


To continue this example we take the ‘single seed’ initial configuration: X =---00100--- and 
look at the time/space evolution of this CA. A calculation by hand shows that this will produce 
the triangle of Pascal modulo 2: 


= 000000 1000000 
F(X) = 000001 0 100000 
F?(X) = 000010 0 010000 
F3(X) = 000101 0 101000 
F4(X) = 001000 0 000100 
F5(X) = 010100 0 001010 
F9(X) = 100010 0 010001 


When computed on a larger scale (Figure 3.2) this reveals a fractal structure also known as 
‘Sierpenski’s triangle’ [59]. © 


3.2 Characteristics of Cellular Automata 


The above example mentioned some terms and characteristics of cellular automata which we will 
now make more explicit. 


discrete space: The space S which determines the dimension and the size of the configurations 
is discrete. In this thesis we will concentrate on the one-dimensional case S = Z (two-way 
infinite) and S = Zp (periodic boundaries with size k). 


finite state set: The state set Q is always finite and typical Q = {0,1} which defines a binary 
CA. 


cell/cell-value: Each space-position i € S identifies a cell which contains a cell-value X; € Q. 


configuration: À configuration X is the concatenation of all cell-values on the space S at time 
t, and therefore an element of Q*?. 


3.3. FORMAL DEFINITION OF CELLULAR AUTOMATA 27 


Figure 3.2: The space/time behavior of the one-dimensional cellular automaton in Example 3.1 
on an initial configuration with only one cell non-zero. The cell-values 1 are colored black, the 
zero-values white; time is going downwards. This self-similar structure is also known as the 
Sierpenski-triangle and has a Hausdorff-Besicovitch or fractal dimension of log(3)/log(2) ~ 1.59. 


local function: The time evolution of the different cells is described by a local function f. Like 
the space of a CA, the time-parameter is also considered discrete. 


neighborhood scheme: The outcome of the local function at position i only depends on a finite 
set of cell-values in the vicinity of i. This set is defined by the neighborhood scheme N of 
the CA. 


global function: Given a space S (compatible with the neighborhood scheme), the local function 
f imposes a global function Fs on the set of configurations. This global function determines 
the space/time behavior of the cellular automaton on an initial configuration. 


The most stringent and typical characteristic of the CA-model is the restriction that the local 
function does not depend on the time t or the place i: a cellular automaton has homogeneous 
space/time behavior. It is for this reason that CA are sometimes referred to as shift-dynamical or 
translation invariant systems. 


3.3 Formal definition of Cellular Automata 


After this informal introduction to one dimensional cellular automata, we are now ready for a 
formal definition of this model. 


DEFINITION 3.1 (One dimensional CA) A one-dimensional, classical cellular automaton is 
defined by the tuple (Q, N, f) with Q the finite state set, N the neighborhood scheme and f the 
local function. The scheme N = {N,,...,N\nj} is a finite set of Z-values with Nj < Nz < -< 
Nin . 


For every appropriate space S this definition gives us a global function Fs : Q? — Q? according 
to: 


Fs(X) = Fs (8x) 


ics 


= QIAN XitNa jX CEN DE) 
ics 
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To simplify the above notation we introduce the following short hands for the terms appearing in 
the local function. 


NOTATION 3.1 If N = (Ni,... , Nigj) is a neighborhood scheme it will be understood that: 
(Xin) = (Kits Xeno Xen) 
and 
(Xa) = (Xa Xa c , Xp) 


This leads us to the the brief and elementary expression: 


"(gx = RFX) 


ies ies 
Notice that a CA F = (Q, N, f) does not determine the space S it is supposed to act on. 


EXAMPLE 3.2 The CA of Example 3.1 is described by the tuple F = ({0,1},(—1,0,1), f) with 
f(,y,2) = @@z. o 


Because the definition of a neighborhood scheme N does not determine the space S it acts on, 
we can apply the same CA F = (Q, N, f) on different structures compatible with the set Z.. Our 
main interest goes to the simple one-dimensional cases Z and Z; which results in the functions 
Fz : Q7 5 Q and Fy : Q* 5 QF with k € Nt. 

If the neighborhood N does not allow any interaction between the cell-values, the CA will have 
trivial behavior. Typically this is the case if |N| = 1. On the other hand it can been shown that 
every CA F can be redefined as a CA with N = {0,1} without loss of generality. This will be 
proven in Lemma 6.5. 


3.3.1 The Cellular Automata Model 


Cellular automata were first used as such by John von Neumann [12, 46]. He showed the possibility 
of a universal two dimensional CA. Although not hard to prove, we will only refer the standard 
literature for a proof that one-dimensional CA can simulate Turing machines and are therefore 
computational universal. Because of their simple and homogeneous and parallel structure, cellular 
automata are frequently used to model physical systems such as gases, liquids et cetera. This, 
in combination with its use as a mathematical abstraction of parallel computation, makes it a 
unique combination of physics and theoretical computer science. The last decade there has been 
a growing interest in theory of cellular automata by computer scientists and physicists because of 
the possibility to simulate large CA with the use of modern computer technology. The articles by 
Stephen Wolfram [58, 59], Tommaso Toffoli [53, 23], Norman Margolus [54], and Howard Gutowitz 
[27] are a good starting point to investigate this modern use of cellular automata. 


3.4 Reversible Cellular Automata 


Given a global function F : Q7 — Q7, we say that the CA F is global reversible if and only if every 
configuration Y € Q? has one predecessor X such that: F(X) = Y. A simple reversible cellular 
automaton (RCA) is described by: k = 2, N = {0,1} and the local function f with: f(xo, 71) = £o 
for every zo,21 € Q. This RCA is just the identity CA with F(X) = X for every X € Q7. By 
this example we see that it is possible to have a RCA with a non-reversible local function (which 
is always the case with k > 2). More complex RCA are also possible but appear to be very rare 
among the set of plain cellular automata. 

The characteristics of reversible cellular automata are extensively described in a review article 
by Tommaso Toffoli and Norman Margolus [27, 52]. For a long time it was unknown if there exists 
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a reversible cellular automata capable of embedding a general purpose computer. In 1976 however, 
Tommaso Toffoli proved that any d-dimensional cellular automaton could be simulated by a d+ 1- 
dimensional reversible cellular automaton [52]. By the existence of a universal Turing machine 
embedded in a 1-dimensional cellular automaton, this proves that there exists a 2-dimensional 
reversible cellular automaton which is capable of doing the same. After that Morita and Harao 
proved the existence of a one dimensional RCA which is universal in its computational power 
(19, 44, 45]. They did this by defining a special subclass of CA which are reversible by definition: 
partitioned cellular automata. This type of CA will also be discussed in this thesis. 

The ‘modern approach’ to cellular automata with the use of computer simulations seems to 
suggest that the old articles on this subject are outdated. Especially for the model of reversible 
automata this not the case. The theory of cellular automata as a computational model is a 
mathematical theory. To avoid ‘re-inventing the wheel’ one should also pay attention to the more 
abstract articles by Serafino Amoroso et al. [2], Gustav Hedlund [28], and Daniel Richardson [48]. 


Chapter 4 


Quantum Cellular Automata 


In this section we will define the model of one dimensional quantum cellular automata. 
This model is a straightforward extension of the classical model that exists for one 
dimensional CA. 


4.1 The Size of the State Space 


Before we describe the various models a note must be made on the size of the QCA we are con- 
sidering. If we take a two way infinite automaton, the number of basis states will be infinite and 
even uncountable. Because we want to describe the general behavior by a transformation matrix 
this is highly impractical. We solve this problem by restricting ourselves to finite, (size k) circular 
bounded automata (also called periodic QCA). This means that the set of basis states corresponds 
with the finite set of functions: Zk — Q instead of the uncountable set of functions Z — Q (Q is 
the set of states for each cell). 

This is somewhat different with the common approach to this problem. Most authors so far 
(John Watrous [57] and Christoph Dürr, Houng Lé Thanh and Miklos Santha [20, 21, 29]) define 
their QCA with a quiescent (non-active) state. Only configurations with a finite number of non- 
quiescent cells are then taken into account. From now on we will only look at one-dimensional 
structures. 


4.2 Quantum Cellular Automata 


With quantum cellular automata (QCA) we refer to the general description of a one dimensional 
quantum cellular automaton. It is a natural extension of the classical definition. 


DEFINITION 4.1 (One dimensional Quantum cellular automaton) A QCA F is defined by 
the tuple: (Q, N, f), where: 


e Q describes the finite set of states 
e N C Z defines the finite neighborhood scheme with N = (No... , Nin} and No < ++: < Nin] 
e f is the local function f : Q > Ho 


In this last definition we see the quantum-mechanical aspect of the automaton. For every k € N+ 
this F gives us a transformation Fy : (3(QF) — €2(Q*) described by 
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cells of QCA at time T 
a a a ae 
cells of QCA at time T+1 
“ae ab alba > 
cells of QCA at time T+2 


Figure 4.1: A typical one dimensional quantum cellular automaton with neighborhood scheme 
N = (—1,0,1). Notice the identical structure when compared to the CA in Figure 3.1. 


for every X € /5(Q*) with ag € C. The values F;|£) are defined as a tensor product of the local 
values £; of the basis state € computed by f: 


Frey = OO Fes mcis E EN) 
JEZE 

= Q Fén) € He)" 
j€Zy 


Note that the addition “j+N” is in Z; and not in Z and also that the range of F(Q") is a subset of 
the proper non-entangled states. The relation between the basis states, the proper non-entangled 
states, the proper states and all possible states in the Hilbertian space is described by: 
k k k 
Q? € (He) G Ho. € &(Q") 


€ = 


Figure 4.1 is a picture of a QCA with neighborhood scheme N = {—1,0, 1}. 


4.2.1 Normalized QCA 


The neighborhood set N defines the possible interaction between the cell values at at each com- 
putational step. Automata with large neighborhood sets allow more ‘internal communication’ 
than QCA with small ‘neighborhood sets’. Before proving that a neighborhood scheme of {0,1} is 
sufficient to implement any possible QCA, we will first have to give the following definitions. 


DEFINITION 4.2 (Contiguous neighborhood scheme) A QCA is contiguous if its neighbor- 
hood scheme is described by a closed interval on Z: N = (No, No +1,... , Nini} 


DEFINITION 4.3 (Normalized neighborhood scheme) A QCA is normalized if its neighbor- 
hood scheme is a contiguous with No = 0, therefore: N = {0,1,... ,|N| — 1]. 


DEFINITION 4.4 (Shift) A shift transformation S on Hox is a transformation, defined by: 


SX) = S| > lal} = So a-Si = o | Qim 


£cQ* £cQ* £cQ* JEL 
for every X € Hox. 


We will often write S|xo,... Z4 1) = |*i,... ,£k-1, £0) to describe the behavior of the shift- 
transformation. By S7 we mean the transformation S, d-times repeated. A negative d indicates 
the use of the inverse translation S1. 
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DEFINITION 4.5 (Shift equivalence of states) Two configurations X and Y € Ho» are shift- 
equivalent if there exists a d € Z such that: X = T*(Y). 


Because we allow negative shifts this is an equivalence relation. 


DEFINITION 4.6 (Shift equivalence of automata) Two QCA F and F’ (both with state-set 
Q) are shift-equivalent if there exists a d € Z such that for every X € Hox we have: F(X) = 
S*(FL(X)). A shorthand for this relation is F = S* o F' or F = F”. 


Note that the value d does not depend on the input X or even the size k of the structure. With 
these definitions we can describe the following lemma. 


LEMMA 4.1 For every QCA F there a exists a normalized QCA F’ which is shift-equivalent with 
F. 

PROOF. (by construction) Given the QCA F = (Q, N, f) we define F” with the tuple (Q’, N', f^) in 
the following way: an equal state set, Q' = Q, a neighborhood scheme N’ = {0,1,... , Nin; - Ni} 


and the local function: f’: QN — He with: 


f' (Xo, Xi... Xii) = f(Xo, XN- Xng—nyy «++ Xp) 


for every X € Q^' With this construction we have F = S o F”. 


This F” is called the normalized version of F. If F' = S40 F we have for every t € N: Ft = Sto F 
(this is because F o S = So F for every QCA F). 


4.2.2 Well-formed QCA 


We start the discussion of well-formed QCA with a counter example. 


EXAMPLE 4.1 Take a QCA defined by ((0, 1), {0,1}, f), with f defined by: 


£(0,0) = |0} (0,1) = 3V2{|0) + |] 
(1,0) = $V2(l0) — |1)} fQ, 1) = |1) 


This gives us the following evolution from the initial state |0, 0,1) € Hos: 
|0,0,1) —r f(0,0) @ f(0,1) & f(1,0) 
z oh & (2o) un } o [58th - my 
= 5 {|0,0,0) — |0,0,1) + |0,1,0) — |0, 1, 1)} 


1 


— rp, 
This last vector is not proper: its norm equals 22/4 ~ 1.17. Because this contradicts the laws 
of quantum mechanics it is said that this QCA is not well-formed: the transformation F3 is not 
unitary. © 


This shows that our definition of QCA allows automata which are not proper. What we mean by 
well-formed QCA is defined as follows. 


DEFINITION 4.7 (Well-formed QCA) A Qca F is well-formed if and only if for every k € Nt 
the corresponding Fy is a unitary transformation. 
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This means that F;, applied to the basis states Q^ has to be both norm and angle preserving in 
the Hilbertian space /3(Q^). Because the range of F;,(Q*) consists of proper non-entangled states, 
it holds for every x and v € Q^ that F,|€) and F,,|v) can be written as a tensor product: 


Fil) = Q zi and F,\v) = Q Yi 
iEZk iEZk 


with xi, yi € Ha. The inner-product of these two vectors can therefore be calculated by the 
multiplication of the inner products of the individual cell-values; 


(Felx), Fk) = (s.s) = [I @.u) 
iEZk iEZk iEZk 


For every x; € Ho we know that (z;,2;) = 1. This means that for every QCA the Fẹ by definition 
preserves the norm of the basis states: 


Uk). Felx)) = (e Ti, @x) = Ius = 1 
i€Zy  i€Zp i€Z 
If we want Fp to be angle-preserving we must have for every € 4 v € Q* (and therefore (J£), v) = 0) 
221999 210) = (e Ti, Qu) = II (sy) = 0 
iCZ,  i€Zp icZy 


This means that there must exist a i € Zp for which (z;,y;) = 0. The following lemma is easy to 
prove. 


LEMMA 4.2 For every QCA F, if F' is the normalized QCA with F = S"(F') (as described in 
lemma 4.1) the well-formedness of F equals the well-formedness of F”. 
PROOF. S is a norm and angle-preserving transformation on Hox. 


4.3 Proving Well-formedness 


In this chapter we will show that the well-formedness restriction can be translated into a local 
constraint. This enables us to formulate an algorithm that decides if a given QCA is well-formed 
or not. 


DEFINITION 4.8 (General function of a QCA) Given a QCA F = (Q, N, f) we define the gen- 
eral function Fz : Q* > (Ho)? by: 


F(X) = &s(Xitw) € (He)” 
i€Z 
for every X € Q*. 


This function describes only the first time step of a QCA on the two-way infinite basis states. Using 
a well-formedness definition for Fz we can investigate the unitarity of the periodic functions Pj 
for all values of k. 


DEFINITION 4.9 (Well-formed general function) The function Fz is well-formed if and only 
if for every X Z Y € Q? there exists an i € Z such that: 


(F2(X)i, F2(Y)i) = 0 
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We will now prove the equivalence of both definitions of well-formedness for QCA. 


LEMMA 4.3 A QCA F is well-formed if and only if its corresponding general function Fz is well- 
formed. 

Pnoor. (If) If F = (Q, N, f) is not well-formed then there exists a k such that Fy is not angle- 
preserving. This means that there exists a x and v € Q^ with y L v and Fy|x) £ Flv). If we 
define the two-way infinite states X, Y € Q7 by: 


X= Q x and Y= Q U 
iez iez 
it follows that for every i € Z we have: 
Fz(X)i = (Felix) moar and Fz(Y)i = (Felv))imoar 


Because X # Y and Fz(X); £ Fz(Y) for every i € Z, X and Y violate the well-formedness 
definition of Fz. 

(Only if) Without loss of generality we assume F to be normalized. If Fz is not well-formed then 
there exist two basis states X,Y € Q7 with Xo Æ Yo and 


(F2(X)i, F2(Y)i) # 0 


for every i € Z. In Appendix B it is shown that with the states X and Y we can construct a 
k < 2|Q[?!N! + | N| and two periodic states X’,Y’ € Q} with |X’) L |Y’) and Fy|X’) £ Fly’. 
This proves that the QCA F is not well-formed. 


In the above lemma, the proof is more important the conclusion: with little extra effort, we can 
deduce the following local constraint for well-formed QCA. This lemma will be the backbone of 
this thesis. 


LEMMA 4.4 For every well-formed QCA F = (Q, N, f) there exist two values p,q € Z with: 


SQFT MD op Sg QPOs Nigam 


such that for every X,Y € Q? with Xo 4 Yo we have: 


q 

[FY = 0 

i—p 
PROOF. Because Lemma 4.3 is ‘if and only if’, it follows (see Appendix B) that for every 
normalized, well-formed QCA F = (Q,(0,...,|N| — 1}, f) there exists a p > -|gp and a 
q < |QU"! + |N], such that for every X, Y € Q7 (with Xo 4 Yo), we have: 


[pX FY) = o 


i—p 


'The constructions in Lemma 4.1 and Lemma 4.2 generalizes this to the desired result. 


One of the important consequences of this lemma can now be proven directly. 


LEMMA 4.5 The well-formedness property of a QCA F — (Q, N, f) is decidable. 
PROOF. If F is not well-formed, there exists a k € Nt such that the finite dimensional transfor- 
mation F; is not unitary. Because k is bounded by 


k < 2|Q| Nin N:+) 


this is decidable in finite time. 
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The importance of of this bounded restriction lies in the fact that it does not depend on the 
local function f of a QCA. For classical reversible CA the range [p, q] is also known as the inverse 
neighborhood of a CA. Previous work by Jarkko Kari [31] suggests that the bound in Lemma 4.4 
can be made smaller. 

If the domain of a well-formed QCA is a proper state-space, the range will also be a proper 
state-space. In that case we can write Fk : Hor — t. If we want to respect the physical laws, 
only a fraction of the possible QCA may be considered meaningful. This resembles the situation 
with classical CA where reversibility is a strong restriction on the automata [2, 48] From now on 
it will be understood that by QCA we mean well-formed QCA. 


4.3.1 A definition of Balancedness 


Recently, Christoph Dürr et al. [21] raised the question about a definition of ‘balancedness’ in the 
case of quantum CA. Here we give a generalization of the classical definition used by Amoroso and 
Patt [2] (which differs from the one used by Maruoka and Kimura [38]). 


DEFINITION 4.10 (Balanced QCA) A QCA F = (Q, f, N) will be called balanced if and only if 


5 Iq, fa? = igit 


rzcQN 
for every q € Ha, 


The following lemma shows the validity of this definition. 


LEMMA 4.6 Every well-formed QCA F = (Q, f, N) is balanced. 
PROOF. If we take k = N\y; — Ni +1, the summation of f(x) is ‘encapsulated’ in the summation 
on the set of basis states Q*: 


316»? = IQ" M? | (a, (5X) ? 


rzcQN XcQ* 


With Lemma 1.1 we can derive the equality: 


QNE S^ He LO) = (aS | $5 eY, cub» 


XEQ*E XeQ* \YeQr-1 


ig" 5 S^ | Y^ |g @Y, (Fax)? 


YeQ*-! VXcQ* 


Because Fj is a unitary transformation such that Hp, (Qr) = Hee with q8 Y a vector with norm 
1 in Hor we reach: 


jos NS] POU CERRO I n oem get oo gem qoe 


yceQ*-! iJ XeQ* yeQ*-1 


which proves the lemma. 


'This lemma tells us that a QCA can only be well-formed if the output values of the local function f 
are uniformly distributed on the state space Hg. This is not a sufficient restriction for properness, 
there exist balanced QCA which are not well-formed (see Example 4.1). 
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4.4 Quiescent Quantum Cellular Automata 


In order to relate this thesis to some earlier results about QCA, we will give a short description of 
QCA with quiescent states: Quiescent QCA. 


DEFINITION 4.11 (Quiescent QCA) A QQCA F is a QCA defined by the tuple (Q, 0, N, f), with 
a quiescent state O € Q, such that f (LIN) = |D). 


This definition is used to overcome the problems of describing the behavior of a QCA in an un- 
countable infinite state space Hgz. The only configurations that are allowed for a QQCA are 
those with only a finite number of non-quiescent states. Because the left and right tail of such a 
configuration (--- ToL, ++ LN ---) will remain quiescent under the action of a QCA, only a 
countable subset of basis states will be necessary to describe the evolution of the system. The set 
of basis states of these finite configurations is denoted by Q*. 

For every QQCA F the time operator on the superposition of finite configurations F, : 02(Q*) > 
£5(Q*) is defined by 


F(X) = F| Y «dà 
EEQ* 


- J aR 


SEQ* 


for every X € L(Q*). The function F.|£) on the basis states Q* is determined by: 


FJS = XO Fian) (€ &(Q*)) 


jEZ 


The well-formedness issue of the function F, is subtle and has some potential pitfalls. Here is a 
first attempt of definition. 


DEFINITION 4.12 (Well-formed QQCA) Given a QQCA F = (Q,0,N,f), the function F, is 
well-formed if and only if for every X,Y € Q*: the norm of F,|X} equals 1 and X L Y implies 
F,|X) L F,|Y). 


Because the domain of the local function f is a subset of Hg the norm of F,(X) will always be 1. 
The following example will show that there exist well-formed QQCA which are not unitary. 


EXAMPLE 4.2 Take the QQCA F = ({0, e), D), {0,1}, f), with: 


f0,0) = f(e,eẹ)=|O) and f(O,e) = f(e, 0) = |e) 


The function F, is well-formed but not injective and therefore not unitary. This is shown by the 
equation F(X) =|... e. ...) which can not be satisfied for X € Ho». o 


'This example indicates that we have to restrict definition 4.12 to the injective QQCA. 


DEFINITION 4.13 (Unitary QQCA) Given a QQCA F = (Q,LI, N, f), the function F., is unitary 
if and only if the function F, is well-formed and for every Y € Q* there exists a X € Hg» with 


Because QQCA are a subset of QCA as defined in Definition 4.1 we can relate the above definition 
of unitary QQCA to that of well-formed QCA. 


LEMMA 4.7 If a QQCA F = (Q,0,N, f) is well-formed as described in Definition 4.7 then the 
F, is a unitary function. 
PROOF. (by contradiction) If the function F. is not unitary there are two possibilities: 
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1. there exist X L Y € Q* with F,[X) L F.|Y) 


2. there exists an X € Q* with no predecessor in the domain Ho». 


Because X,Y have a left and right tail which is quiescent, it is possible to take k € N* sufficiently 
large, such that the same possibility holds for Fk. This shows that Fẹ is not unitary and therefore 
proves F not to be well-formed. 


The converse of this lemma would be: “If a QQCA has a unitary F}, then the QCA F is well-formed.” 
The following example shows that this is not the case. 


EXAMPLE 4.3 Take the QQCA F = ((L1,0, 15, L1, {0,1}, f), with f defined by: 


f(,0) = [p f0, x) = |D) 

f(z,D) = |x) f(z,y) =| e y) 
for every x,y € {0,1}. Although the function F, is unitary, the QCA F = ((D,0, 1), N, f) is not 
well-formed. This is shown by F5|0,0) = Fə|1, 1) = |0,0). © 


We can summarize the above lemma’s and examples in the following hierarchy with F a QQCA: 


Fk is unitary for every k € Nt 
Fz is well-formed 


F, is unitar F, is well-formed 
f= y f= 


It is still unknown how to define well-formedness for non-quiescent state QCA on a two-way infinite 
structure with an uncountable dimensional state space Hgz. 


4.5 Partitioned Quantum Cellular Automata 


The results of the previous sections showed how to decide if a QCA is proper or not. What we 
still do not know is how to construct proper QCA with non-trivial behavior. This situation is 
no different with the problem of non-trivial reversible classical CA. A solution to this problem is 
(partly) given by partitioned CA [19, 44, 45]. Partitioned CA are constructed in such a way that 
the well-formedness of the automaton is implied by their definition. John Watrous applied this 
method to construct a QCA which can simulate a quantum Turing machines [57]. 


DEFINITION 4.14 (Partitioned QCA) A partitioned quantum cellular automaton (PQCA) is 
defined by a tuple (Qr, Qc, Qn. g), with Qr, Qc and Qn finite state sets and g : Ha,xQcoxQrn — 
HQLxQcxQr 8 unitary transformation. 


A PQCA defines a special type of QCA F = (Q, (-1,0, 1), f) with Q = Qr xQcxQr. This product 
state set reflects the idea to decompose the cell-values into three disjoint parts Left, Center and 
Right. Each state x € Q is therefore identified by three sub-states: x = x! x z^ x x". This is used 
to translate the unitary function g into the local function f : Q? — Hg by: 
f(z-1,20,21) = g(a) x 25x x”) 

for every z.1,29,2, € Q. Because g is a unitary transformation, the QCA F will have a local 
behavior which is also unitary. From this it follows that F' will be a well-formed QCA. See Figure 
4.2 for an example of such a PQCA. With the use of PQCA we can define QCA in a straightforward 
way without having to worry about the well-formedness constraint. This has proven to be a fruitful 
concept in the theory of reversible CA and proper QCA. It assures us that the model of QCA is a 
non-trivial way of describing parallel quantum computation. 


LEMMA 4.8 There exists a proper QCA which can simulate a Universal quantum Turing machine 
with polynomial time complexity. 

PROOF. See the original paper by John Watrous [57] which combines the earlier results on re- 
versible CA and quantum Turing Machines. 
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Figure 4.2: A partitioned quantum cellular automaton. The local function corresponds to a unitary 
transformation on the three sub-parts of a cell-value. The well-formedness of the local function 
gives us a global behavior of a well-formed QCA. 


Chapter 5 


Quantum Gate Cellular Automata 


The definition of the quantum gate cellular automata (QGCA) is strongly related to 
the notion of quantum gate arrays. By replacing the local function of a QCA by a 
proper quantum gate we get a QCA which is proper by definition. This resembles the 
construction of partitioned QCA. Because of its central role in this thesis an extensive 
description will be given. 


5.1 Description of QGCA 


A quantum gate cellular automaton (QGCA) is a parallel and uniform circuit with identical quan- 
tum gates. See Figure 5.1 for an example of such a structure. 

To avoid the problems of infinite sized systems, we will only consider QGCA with a circular 
bounded structure. The number of input/output values will be called the fan of the gates. For a 
proper quantum gate the fan-in has to equal to the fan-out, so no confusion is possible with this 
terminology. If we have k gates with fan n the number of cell-values will be nk. For clarity we 
use two indices to indicate an individual cell on the structure Zp. Every basis state € of such a 
system can be described by the tensor product of the individual cell values: 


E) = 28 @2,@--- O20 O21 @-:-@29@---@at] 


which equals 


e - @(@+) - Qz 
j€Z& \i€Zn JEZR 
with a^ € Q for every i € Zn and j € Zy. 

Because the gates are in a parallel position, the behavior of a layer of gates is described by 
a unitary matrix U which is the k-fold tensor product of the individual gates M. Therefore 
U = MW. 

The neighborhood scheme describes the ‘wiring’ between the gates. Each gate has the same 
neighborhood scheme which does not depend on the size k of the circuit. 


DEFINITION 5.1 (Neighborhood scheme) A neighborhood scheme P for gates with fan n € 
Nt is defined by a tuple (n,o, 9). Both o : Zn — Zn and $ : Zn — Z are functions, with o a 
bijection on Zp. 
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Figure 5.1: A part of a quantum gate cellular automaton. All the gates M are identical. The 
neighborhood scheme P defines the wiring between the gates. Because M is proper gate, the 
global transformation function of this system will also be proper. We therefore know that this 
structure resembles a well-formed QCA. 


A neighborhood scheme P permutes the values of a basis state £ € Q"^ in the following way. 


rie = P(@ (Bled) 


j€Zyk \i€EZn 
®&) (& 9) 


jeZy Ni€Z, 
This shows that P corresponds to a function Znk — Z4; and also that the description of P does 
not depend on the involved set of states Q. Because o is a bijection, P will also be a bijection for 
every k € Nt. This is necessary because we want P to be a proper transformation. If P acts on 
a superposition, the behavior of P will respect this superposition: 


PX) = P| Yo ad] = So aP 


£eQ"nk £eQ"nk 
for every |X) € Hane. 


EXAMPLE 5.1 The neighborhood scheme as shown in Figure 5.1 is described by (n = 3,0, ¢) 


with: 
c(0 = 2 $(0 = -1 
oll) = 1 and 4$(1 = 0 
c(2 = 0 (2) = 1 


© 


The reader is encouraged to check the above example, to get a better insight in the used formalism. 
Note that the neighborhood scheme acts before the layer of gates is applied. 


5.2. PRELIMINARIES Al 


5.1.1 Formal Definition of QGCA 


We are now ready to give a formal definition of a QGCA. 


DEFINITION 5.2 (Quantum gate cellular automata) A QGCA F is defined by (Q,n, M, P) 
with Q a finite state set, n € Nt, M : Han — Han a proper quantum gate with fan n and P a 
neighborhood scheme (n, oc, ¢). 


As with QCA this definition does not specify the size k of the array. A layer of a QGCA corresponds 
to the sequential combination of a neighborhood scheme P and a row of gates Ml. One layer 
can be identified by one time step of the QGCA. The global behavior of a QGCA F — (Q,n, M, P) 
of size k is denoted by F and obeys: 


FX) = MUoP| Y. a 
&cQn* 
= M aM o PIE 


£eQ"nk 


for every X € Hans. The effect on the basis states € is defined by: 


M™ oP = M'"oP| R (8s) 


jEZk \tEZn 


II 


Le (eet) 


j€Z& Ni€Z, 
- eG) 
- &(v(&«t.)) 
j€Zk i€Z, 


The function of a sequence of t layers Fẹ will be described by F¢. 


5.2 Preliminaries 


Before continuing with the general theory of quantum gate arrays, we will take a closer look at a 
particular neighborhood scheme which plays an important role: the Shift neighborhood scheme. 


5.2.1 The Shift Neighborhood Scheme 


In this thesis we will concentrate on the Shift neighborhood scheme S, which simply ‘moves’ all 
the values one cell leftwards (see Figure 5.2 for an example). With the fan-value n this gives us 
the tuple (n,i + 1,15 1(2)) with: “i +1” an addition in Z, and 1, the function with 1e(x) = 1 if 
x = c and 1,(x) = 0 otherwise. 


EXAMPLE 5.2 Because of their size, it is not practical to actually write out the transformation- 
matrices. As an example of the matrix S for the simple two state system with kn = 3 the equation 
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b Se y See ---- MM H kn cell values at time T 
"A [ y^ Shift-neighborhood-scheme 


gate #1 gate #2 -—— gate #k k gates circular bounded 
ERES [| kn cell values at time T+1 
2 " P 4 S S Shift-neighborhood-scheme 
Y Y 
gate #1 gate #2 mma gate #k k gates circular bounded 
y ! 279 UD ND VV 


Figure 5.2: The shift neighborhood scheme. The neighborhood scheme S describes a leftwards 
translation on the cell values. 


S|0, 1, 1) = |1, 1, 0) looks like: 


10000000 ones Ojo00) 

0-0 0 0 1 205.00 Biby Ojoo1) 
01000000 ices) Ojor0) 
|00000100 1011) — | 90m» å 
ee ee ee ee A Geach Fe retos | rien 
00000010 Gites Oj101) 
00010000 (iftos liio) 
00000001 Öar Otai 

© 


From now on when using S in the context of QGCA, we will mean the shift neighborhood scheme. 
The inverse of S (a shift rightwards) is denoted by S^. 


5.2.2 Two state QGCA 
The following lemma shows us how to simulate a QGCA with an arbitrary state set Q by a QGCA 
which uses qubits. 
LEMMA 5.1 Every QGCA F — (Q,n, M, P) can be simulated by a two state QGCA 
G E (10, 1}, n’, M’, P’) 


with an equal number of computational steps. 

PROOF. (by construction) First expand the set Q to a set R D Q such that |R| = 2^ with 
à = [logs |Q|]. Now we can use a bijection ¢: R — {0,1}* which, induced to the n cells of each 
gate, gives us a bijection ® : R” — {0,1}™. With this ® we define the new quantum gate M’ on 
Hyo,1} by: for every x € R” by: 


; " $ !(M(r) ifre Q” 
M'($(r) = eee otherwise 


for every x € R”. The new neighborhood scheme P’ = (nA, o', ¢’) is obtained by splitting each 
‘wire’ into A parts: 


o'(iA +t) = a(i) and  '(ià+t)=¢(i) 
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for every i € Zn and t € Z). 

This QGCA G = ((0,1], nA, M’, (nÀ,c', ¢’)) simulates the original QGCA according to the 
bijection $. If an illegal state is encoded (x d Q”) the M'-gate behaves like the identity. This 
shows that M' is a proper transformation. 


5.2.3 Combining neighborhood schemes 


Here we will show how to express a neighborhood scheme R o P which is a combination of two 
initial schemes P and R. Given the two neighborhood schemes P = (n,o, 9) and R = (n,<, p}, we 
can define the action of P by: 


r(@(@x)) - e(ess.) - @(@x) 
jEZk Ni€Z, j€Zy& \iEZy j€Z& Ni€Z, 


If we use the y expressions to describe the effect of R4 we get: 
i = s(i) 
n(@(@x)) - @(@x%o) 
j€Z, \iEZn j€Z. \i€Zn 


a(i) ): 


The combination of R and P therefore obeys (using the equation y = Tjj) 


RoP| C9 (& J = & (& o 
jeZ. NieZ, j€Z. Ni€Z, 
'This shows us that we can write: 
RoP = (nsi) e)o(no()60) = (mooc(i) p(t) + 6oc(0) 
With regard to the shift neighborhood scheme we have: S = (n, i + 1,14 1(;)) and 


S"oP = (nil)o(no(i)9() = (mo()l-c-e(i)) = Pos” 


5.3 Every QGCA describes a QCA 


If we look at the figures 4.2 and 5.1 we see a great resemblance. We already know that PQCA are 
a special kind of QCA. We will now prove the same holds for QGCA. 


LEMMA 5.2 Every proper QGCA F = (Q,n, M, P) defines a well-formed QCA QCA F = (Q", N, f). 
PROOF. The correspondence between the values of the QGCA and those of the QCA is expressed 
by: 

s = odes 


for every x; € Q”. The neighborhood scheme P = (n,c,$) gives us the neighborhood set of F 
according to: N = {@(i) | i € Zn}. The local function f : (Q")" > Han now obeys: 


fey) = M (e«t) 


iEZn 
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By this construction, for every k € Nt it holds that: 


«[&(&«)) = »[e(e«) 


jEZk \tE€Zn j€Zy& \t€Zn 


= [e (ts) 


jEZk \tE€Zn 


- &(v(&«u.)) 


j€Zx 1€Zn 


= Q) f(x.) 


j€Zk 


= Fk Q Tj 


JEZE 


As with PQCA the well-formedness of a QGCA (Q, n, M, P) is certified by the well-formedness of 
the gate and the neighborhood scheme. 


Chapter 6 


QCA, PQCA, and QGCA are 
Equivalent 


In this chapter will prove the equivalence of all three models for quantum cellular 
automata. By ‘equivalence’ we mean that every QCA can be simulated by a QGCA, 
every QGCA by a PQCA, every PQCA by a QCA etc. All these simulations have at 
most linear slowdown. Our main problem is to prove that every proper QCA can be 
simulated by a QGCA. 


6.1 A Definition of Equivalence 


Suppose we have two QCA F and G. If we say that G simulates F, we mean that there are two 
‘simple transformations’ $, and V, such that ®; o G% o V, equals F*. The G>-expression indicates 
a simulation with only linear (proportional to A) slowdown. The transformations must be ‘simple’ 
because we want G to do the actual calculation. Or in the words of Marvin Minsky [43]: 


( ... ) for if one is permitted an arbitrary partial-recursive computation to do the 
encoding ( ... ), then one could use as the code the result of the Turing-machine 
computation itself, and this would surely be considered a cheat! 


With this warning in mind we pose the following (inductive) definition: 


DEFINITION 6.1 (Simple transformation) Let P,Q, R be state sets and a,b, c, d € Nt. Sim- 
ple transformations depending on t € Z, are constructed by the following three procedures: 


1. If à: P* — QP is a total function, then (t) : P*” — Q” defined by: 
e(t) (& Jj = Q [Orie es ERU 1) lic m er 
iez i€Z 
is a simple transformation. 


2. If ®(t) is a simple transformation, then the transformations S o ®(t), S^! o P(t), St o ó(t) 
and S^! o ®(t) are also simple (with S the shift translation). 


3. If ®(t) : P^? — Q% and Y(t) : Q** — R? are simple transformations, then the combination 
O(t) : PeZ — R described by O(t) = V(t) o ®(t) is also a simple transformation. 


A simple transformation (t): P7 — Q?Z defines for every k € N* the simple transformation 


$, (t) : TL psa — Hore which obeys the superposition principle of the state space. 
Now we can formulate the definition for simulation of quantum cellular automata. 


45 
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DEFINITION 6.2 (Simulating QCA) A QCA G — (P, N', g) can simulates a QCA F — (Q, N, f), 
if there exists a tuple (A, , Y) with e(t) : Q% — P'7 and W(t) : P2 — Q™ simple transforma- 
tions and À € N such that for every t,k € Nt and X € Hora we have: 


Fha(X) = ®ra(t) o Giro Vis(t)(X) 
Or in other words: 
Ffa = Oralt) o (Gr) o Vit) 
If no misunderstanding is possible we can shorten this to: Ft = $, o G o Wy. 


When F can be simulated by G, we will write “{F} < {G}”. This relation is certified by (A, ®, V), 
this tuple will also be called the encoding of F into G. 
With this definition we can define the notion of equivalence for quantum cellular automata. 


DEFINITION 6.3 (Equivalent QCA) Two Qca F and G will be called equivalent if and only if 
F can simulate G and G can simulate F. 


Consequently we will denote this by *(F) ~ {G}”. For every Qca F, G and H we now have: 
1. {F} 3 {F} 
2. If {F} < {G} and {G} < {H}, then {F} < {H} 
3. If {F} < {G} and {G} < {F}, then {F} ~ {G} 


If we take the ‘identity-encoding’ (1,7, I), we see that 1 holds, whereas 3 holds by definition. The 
second rule needs more explanation. 

If {F} < {G} is certified by (A, ®, V) and {G} < {H} by (u, A, ^), we have (in shorthand): 
Ft = 6, o G^ o V, and G) = A; o H!? o A; and therefore: 


Ft = G0A,,0H oho, 


By defining =; = 9, o Ay, and Y, = Ay; o V,, we see by Definition 6.1 that = and Y are simple 
transformations. This shows that (uA, E, Y) certifies the {F} < {H}-relation. 

Because the relations “<” and “~” are defined on sets of QCA we can extend their domain in 
the following way: 


DEFINITION 6.4 (Simulating sets of QCA) If.A and B are sets of quantum cellular automata, 
we say that B can simulate A if and only if for every F € A there exists a G € B which simulates 
F. This will be denoted by “A < B". 


DEFINITION 6.5 (Equivalent sets of QCA) If A and B are sets of quantum cellular automata, 
we that A and B are equivalent if and only if both A < B and B < A holds. This equivalence is 
denoted by A ~ B. 


With this definition we have for every set of quantum cellular automata A, B and C: 
1. AZSA 
2. If A< Band B <C, then A <C 
3. If AX B and B < A, then A~ B 
4. Ac A 
. Az B, then B~ A 


oO 


6. If A~ B and “B ~C”, then AYC 


6.2. SOME PRELIMINARY RESULTS 4T 


These rules show that < is a partial ordering and ~ is a equivalence relation. A few concluding 
remarks about the above definitions: 


e A~B does not necessarily mean that for every F € A there exists a G € B, with F and G 
equivalent. 


e The constant “A” indicates the linear slowdown of the simulation. 
e If ACB then automatically A < B. 
e If F is shift equivalent with G then F and G are also equivalent. 


e An interesting question is whether the <-relation is total or not. Does for every QCA F and 
G it holds that {F} < {G} or {G} ~ {F}? 


6.2 Some Preliminary Results 


The strategy of this chapter is to reduce the different sets of QCA, PQCA and QGCA to each other 
in a step-by-step way. Let us first summarize the relations we already know. We begin by using 
the fact that QGCA and PQCA are subsets of QCA. 


LEMMA 6.1 PQCA < QCA. 
PROOF. Partitioned QCA are just a special type of QCA, therefore PQCA C QCA. 


LEMMA 6.2 QGCA < QCA. 
PROOF. Lemma 5.2 proves QGCA C QCA. 


To reduce the class of QCA we use the following lemma. 


LEMMA 6.3 The subset of normalized QCA is equivalent with the general set of QCA. 
PROOF. Normalized QCA are a subset of QCA. Lemma 4.1 proves that every QCA is shift-equivalent 
with a normalized QCA, thus QCA ~ QCAnormalized: 


Now we will prove a new result which will simplify our reasoning about simulating QCA consider- 
ably. 


LEMMA 6.4 The set of QCA with neighborhood set N = {0,1} is equivalent with the set of 
normalized QCA. 

PROOF. (by construction) For every normalized QCA F defined by (Q, {0,... ,|N| — 1}, f) there 
exists a QCA G with neighborhood set Na = {0,1} which simulates F with no slowdown. In 
order to get a neighborhood size 2 we expand the set of states: G is described by the tuple 
(Qt, (0, 1), g). The local function g : (QIN!-1)? — Hgin|-1 in this definition obeys 


g (TN) TIN [:(2|N|-2)) = fau) 8 Faq) 9 7:9 f(@qnj-1)x21N|-2)) 


with x; € Q and therefore £a:(a+|N|-2) € QU, 
The encoding of F into G is established by the simple transformations: 


y : (Q)INI-1 > (QNI 71) with — W(zi4w|-1) = 21 8: @2nj-1 
$: (QUU yes (Q)NI-1 with (21 Q- Q ziv| i1) = zi4N|-a1 
The tuple (1,6, V) now certifies {F} < {G} with Ng = {0,1}. 


A simple extension of this lemma will be the last of our preliminary results. 


LEMMA 6.5 The set of QCA with neighborhood N = {0,1} is equivalent with the set of QCA. 
PROOF. Combine the results: QCA © QCAnormalized ANd QCAnormalized ^ QCAN={0,1} 
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6.3 Shift Neighborhood Scheme is Universal 


In this section we will prove that every neighborhood scheme can be simulated using the shift- 
scheme S. It will be understood that k is the number of gates, n is the fan of the gates, P and R 
are neighborhood schemes and M, M2,... are proper quantum gates. 


6.3.1 Periodic QGCA 


A periodic QGCA is a non-homogeneous QGCA. It is build by a periodic pattern of different gates 
and neighborhood schemes. 


DEFINITION 6.6 (Periodic QGCA) A periodic QGCA F of p layers is defined by a tuple 
F = (Q, n, u, Mo, ... ,My-1,Po,.-- , Pa—1) 


with M; proper quantum gates and P; proper neighborhood schemes for O < i < u. The unitary 
operator of the first u layers of this QCA is defined by: 


m = I (ana) 


(Notice the order of the index i and the fact that Fẹ is defined as the operator of several layers.) 


EXAMPLE 6.1 With the 3-periodic QGCA F = (Q,n,3, Ao, A1, A», Po, Pi, P»); the behavior of 
the first 3j layers is described by: 


; E 
Fi = (abl. Pa. A}. B AD. Po) 
© 


A periodic QGCA with u layers will be call a p-periodic QGcA. An 1-periodic QGCA is again a 
regular QGCA. This definition is a side-step from the regular classes QCA, PQCA, and QGCA to 
allow a more flexible way of defining quantum cellular automata. After that it will be proven that 
every periodic QGCA corresponds to a QGCA. 


6.3.2 Simulating Neighborhood Schemes 


By allowing different layers in a QGCA we can construct a periodic QGCA that mimics the behavior 
of a neighborhood scheme. This is shown in the following lemma. 


LEMMA 6.6 For every QGCA F = (Q,n, M, P) there exists a periodic QGCA G defined by 
G = (Q, n, u, Mg «s ,Mp-1, 8.0259) 


which is shift equivalent with F. 
PROOF. See Appendix C for a proof by construction of this lemma. 


Figure 6.1 shows us an example of such a simulation by a periodic QGCA. 


EXAMPLE 6.2 The QGCA F = (Q,3, M, P) with a neighborhood scheme P = (3,0,¢) and 


c(0) =2 (0) = —1 
o(1) 21 o(1) 21 
o(2) =0 (2) =0 
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original neighborhood scheme 


simulation by permutation gates 


= GHEE 


Figure 6.1: Simulation of the neighborhood scheme P in Example 6.2. The two additional layers 
of the periodic QGCA G consist of permutation gates which mimic the wiring of the original 
neighborhood scheme. The last permutation gate is combined with the gate M of the QGCA F. 
Because the periodic QGCA only uses the Shift scheme this simulation will will be equivalent with 
initial QGCA F up to a S? translation for every layer of F. 


can be simulated by a 3-periodic QGCA 
G = (Q,3,3, Mo, Mi, M5, S, S, S). 
The gates of this automaton are described by: 
Molz, y, z) = |x, y, 2) Mi|z, y, z) =|2,2,y) Map|, y, z) = M|x, z, y) 


for every z,y,z € Q. See Figure 6.1 for an illustration of this simulation by G. Notice how the 
last permutation gate and the M gate are combined to one gate M». o 


6.3.8 Simulating Periodic QGCA 


The use of periodic QGCA will only be temporary because we can simulate any periodic QGCA by 
a regular QGCA. 


LEMMA 6.7 For every periodic QGCA F = (Q,n, u, Mo, ... , Mu 1, $,... , S) there exists a QGCA 
G = (Q,n + A, M', S) with A = [log(u)/log(|Q|)] which can simulate F. 

PROOF. Let '..." be an injective function: Z, — Q^. Define M' : Han+a — Hanta for every 
j € Zp by: 


M86, vic itat j"; £n-1) = MG esu ,£n-1) 8| +10) 
We will use the simple transformations: V : Q” — Q"^^ and ®: Qrt+4 — Q” with 


V(zo,... xac) = (zo,... a1) & ("0") 
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and 
$(rg,... Zu 1,4 1,... d A) = — (x9,... tni) 
Now (A = u, ®, V) certifies the {F} < {G} relation according to: 


OA) e Giint A) Q Ven = Fin 


An example will show us how this procedure can be applied to a small periodic QGCA. 


EXAMPLE 6.3 If a 2-periodic QGCA F = ({0,1},3,2,Mo, M1, S, S) has the following behavior 
for the circuit size k = 2: 


F(X) = M? . S- MË! - S|a0, x}, £2, 0°, 21, 22) 
MIL CMT E d) 

= My" -S (Mold, 25,21) & Moli, 21. 20)) 
= MÝ -S (|y. yo: vo: Vb Vo Yi) 

= MP (yj vo 9o; Vi Y? 9d) 

= Milyo yo Y?) D Milyi vi Yo) 

= |2020 20: Z1 25 21) 


= (2) 


Then, with 
U(xp, To T0) = (20 2020; 0") and B(x}, 29,29; T'o) = (20 20 0) 


the tuple (A = 2, 6, V) certifies the simulation of F by G = ((0, 1), 4, M', S) as is illustrated by: 


2 
Ps o (wu? 8) o W6(X) = Bg o M’? S. M'Pl - S|x§, £929. "07, £f, £1, £7,707) 
= ds o M? S- M'P |æ}, £2,707, 29,23, 22,707, 22) 
= $a o M"? - S (M'|x$, 29, 0, 21) & M'|x1, x1, 0, 25)) 


= $5oM'" .S(Mo|zj, 20,22) & | 17) & Molz}, x7, 28) € | 1) 
= GoM Syd yl, yo, 1' y ul y 1) 

= go M'" ly yo, Ty yl yo 1D) 

= $5oM"yj,yo, 1 91) D M'|yi, yi, 1. yo) 

= 430 Mi|[yo; yo V1) Q | 07) & Mi lyi, vt vo) @ | 07) 

— 6s|29, 29,20 07, 27, z1, 21, 0 ) 

= [20,201 20:20 20 71) 


- (Z) 


o 


With the use of above lemmas we can now state the following theorem about the simulation of 
periodic QGCA. 


THEOREM 6.1 Every periodic QGCA F = (Q,n, p, Mo,... , My 1, Po, ... , P, 1) can be simu- 
lated by a regular QGCA G = (Q,n', M', S) with only linear slowdown. This can be abridged by 
the statement: QGCA periodic ~ QGCAg. 
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PRoor. If apply the proof of lemma 6.6 for every neighborhood scheme P; we have a periodic 
QGCA which only uses the shift neighborhood scheme. With lemma 6.7 we can construct a regular 
QGCA G which has the same behavior and uses the shift neighborhood scheme. 


This is a strong equivalence relation within the class of QGCA. We are now able to define our 
QGCA in terms of periodic QGCA because we know that eventually we can change those into a 
regular QGCA again. 


6.4 Every QCA can be simulated by a QGCA 


Using the above results we will prove that every well-formed QCA can be simulated by a QGCA 
with only linear slowdown. 


LEMMA 6.8 Every proper QCA F can be simulated by a periodic QGCA G. 
PRoor. We already know by Lemma 6.5 that we can restrict this proof to the QCA with neigh- 
borhood set N = {0,1}, we therefore assume F = (Q, {0, 1}, f). 

Lemma 4.4 shows us that there exist a < 0 < b € Zand a function g : R > Q with R C i-a 
which calculates the local inverse of F. This function obeys: € L x if g(£) Z g(x) for every £, x € R. 
The function g is called the local inverse because with x; € Q and |y;) = f (xi, £i+1) the function 


lUa; Vai; --- ;y-1; T0, Y0, Urs --- Yo) Yu 8 lUa; Vac; -- >Y-1;Y0, Y0; Y1,--- 57») 


defines a proper transformation which respects the inner product. This proves the well-formedness 
of the transformation 


|a sus :J—1; T0; Uo Yl; , b) — fo [as sss : Y0; Y0; Y1,- , Yb) 
= [Ya .. ,y-1) 9 {a0, 0) + 8/1, @ [yr ... , yo) 


for every x; € Q, and f(x;, 2441) = |yi) with f(x, 71) = |yo) = a|0) + 8|1). 

We define the following A, B and C gates and their ‘requested behavior’ (see Lemma 2.1) with 
£j € Q; yj = f(zj, cj41) € Ha and r 2 6 — a-- 1. The A-gate calculates the r — 1 values of y on 
the left-side: 


r—1 r—2 
Al Qirra] = | lesa) | @ ltr- £1) 
j=0 j=0 
The B-gate calculates the remaining value of y: 
r—2 T 
B||69lz.v) | @ler-1,2r-1) 8919) |] = Riley) 
j=l j=l 


We use the C-gates to replace the z values by the new y values. Each C gate replaces one x value. 
The well-formedness of this transformation is affirmed by the existence of the inverse function g. 
For every 0 € t < r — 1 and z; € (z;, yi] and |yi-4) = o0) + 81): 


t-a-1 t+r-1 
C, Q [Zj Vj) & eas Yt-a) & Q 123, yj) 
j-t j=t—atl 
t-a-1 t+r—-1 
= 69 [Zj Yj) | 9 luta; Yt-a) ® 69 [Zj uj) 
j-t j=t—atl 


t+r-1 


i—a—1 
Q lz; yj) | & {al0,0) + Bll, 1)} 8 Q 125,45) 
j=t 


j—t—a41 
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By the indices of the x and y values in the above definition it follows that: 


r—1 
II (CF - 82) s?e-9. p . 82. AM. S| Q Inu) 
t=0 


j€Zkr 
r—1 
k 
= II (ci | Js?) Q |Er2+j, Yr—2+5) 
t=0 j€Zk. 


= © luv) 


JEL er 


= Q |f (25, mici), f (£i, £i+1)) 


JEL Kr 


By Lemma 2.1 and the existence of the function g we know that A, B and C are all proper gates 
(note that all C-gates are identical). The proper (r + 2)-periodic QGCA 


G = (Or 2 MB ee , Co, S2, 82, 92079. 95:792) 


therefore simulates the QCA F. This is certified by the encoding (1, ® : Q” — Q", V : Qr — Q?n) 
With (vos odds) (£0, £0, £1, Viger. ,Zr—1;, &r—1) and O=— V4, 


We use this lemma to proof that the class of QCA can be simulated by the the class of QGCA. 


THEOREM 6.2 Every well-formed QCA F can be simulated by a QGCA G. 
PROOF. A well-formed QCA F can be simulated by a periodic QGCA (see the above lemma) and 
every such periodic QGCA can be simulated by a QGCA G (Theorem 6.1). We therefore know that 
there exists a QGCA G with {F} < {G}. 


6.5 QGCA and PQCA are Equivalent 


Because every PQCA is à QCA (which can be simulated by a QGCA), we only have to prove that 
every QGCA can be converted into a PQCA with the same behavior. Again we will restrict ourselves 
to the shift-scheme which covers all possible QGCA. 


LEMMA 6.9 Every proper QGCA F can be simulated by a PQCA G. 
PROOF. The QGCA F is defined by (Q, n, M, S). We construct a PQCA G described by 


G = (Q' = Q, x Qc x Q5. 3, {-1,0, 1}, 9) 
with Q; = Q, Qo = Q"^!, and Qk = {nill}. The simple transformations are V : Q” — Q 
and ® : Q' — Q”, with V(xo,... ,z4 1) = Vo X z1a-1 X nill and $ = Y-t. The local function 
g : (Q)? — Ha obeys: 
g(z,y,2) = g'(zr xyoxzg) = M(yoxzr)xnil € HoxqQn-1x(niilly 


with x,y,z € Q'. Because M is a proper gate, g’ will be a unitary transformation. G is therefore a 
proper PQCA, which simulates the original QGCA with no slowdown. This is certified by (1, ®, V), 
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as is shown by (t — 1 and k € N*): 


P o Gk o Upkn Q (8) = Dr o Gr Q) (aj x æ} xx 277 x nill) 


jEZk NiCZn jEZk 


DO) E eta) somni) 


j€Zx 
Bo 1 n—1 „0 
= Q) M(zj,... 2% ; 541) 
j€Zx 
- ws[e(e«) 
jEZk \t€Zn 


Because QGCA ~ QGCAg we may conclude: QGCA ~ PQCA 


6.6 Conclusion 


THEOREM 6.3 QCA ~ PQCA ~ QGCA. 
PROOF. A summary of the results of this chapter gives us: 


1. QCA < QGCAperiodic (Lemma 6.8) 
2. QGCAperiodic X QGCAg (Theorem 6.1) 
3. QGCAs < PQCA (Lemma 6.9) 


4. PQCA < QCA and QGCA < QCA (Lemma 6.1 and Lemma 6.2) 


The theory of quantum gates and quantum networks assumes: Q = {0,1}. For this reason we 
state the following additional lemma: 


LEMMA 6.10 Every well-formed QCA F = (Q,r, N, f) can be simulated by a QGCA G = ((0,1], n, M, S). 
PROOF. For every QCA F there exists a QGCA F’ = (Q,n', M', S) with {F} < (F'). F’ can be 
simulated by a QacA F” = ((0, 1}, n", M", P) (lemma 5.1). By theorem 6.1 there exists a QGCA 
G = ((0,1), n, M, S) with (F") < {G}. Therefore: {F} < {G}. 


This last lemma enables us to use the known results on quantum gates [4, 16, 18, 60] when 
analyzing QCA, without loss of generality. By doing this, the next section shows us the existence 
of a universal QCA and answers a question about simulating QCA on Quantum Turing Machines 
[21]. 

The above theorems and lemmas also hold when restricted to classical CA. "Theorem 6.3 
therefore proves (in the case of 1d-CA) the conjecture made by Toffoli and Margolus [55] about 
the structural invertibility of (classical) reversible CA. An independent proof of this conjecture 
has been made by Kari [32], which also holds for 2d-CA. 


Chapter 7 


A Universal Quantum Cellular 
Automaton 


In this last chapter we will prove that there exists a QCA which can simulate any 
other QCA with only a linear time slowdown. This QCA will therefore be called a 
Universal quantum cellular automaton. 


7.1 Simulating QCA with QTMs 


With the results of the previous chapter, we are now able to solve a well-known problem: “Is it 
always possible to simulate a QCA on a quantum Turing machine?” The answer will be affirmative. 


LEMMA 7.1 For every well formed QGCA F = ({0,1},n, M, S) there exists a well formed QTM 
T = (£, K,ô) which simulates F. The simulation of the function FL will have O(||) space com- 
plexity and O(||LI) time complexity. 
PROOF. Let the alphabet set © equal (0, 1,0}. Given the unitary transformation M on 7149,13», 
we first define two QTMs: Ts and Ty. 

The Ts simulates the shift neighborhood scheme of the QGCA, it is defined by: 


r z 0 1 n—1 s ia zb 1 2 n—1 0 
|" start"; O, £0, £o- 2g ,,Lh:::) — mg | halt"; O, £0, £0- s Eki to Ot) 


The time/space complexity of this well formed QTM will be O(k). The description of Ts does not 
depend on k (we want T to be valid for any k). 

The second QTM, Tm, has to simulate the unitary transformation M'*! which corresponds to a 
concatenation of k transformations M. Its behavior is described by the following sequence which 
will be explained afterwards. 


[k] 


| start |; Mor YD o ls So) on, ['start!; ) & M|yo) & |y], 5 ec) 
ry |'start'5 D) e Mlyo) e Mly) 9 [yo D) 
Ty, “start '; ) ® Ml|yo) -: & M |yx 1) 8 |, -- ) 
— my, |'halt"; O) & Mlyo) --- & M |yx-i) @ |, --+) 


This QTM applies M to the first n bits of the input string after which the head moves n places 
to the right. This process is repeated until the L-symbol is read. Now the head returns to the 
leftmost L)-symbol and the QTM halts. Because M is a unitary transformation such a well formed 
Tm exists (Lemma 2.3) Again the QTM is independent of k and has time/space complexity of 
O(k). 
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This shows that Tm o Ts simulates one layer of the QGCA F. By repeating this algorithm 
t times we simulate the global function F¢ of the QCA F. Because the time complexity will be 
O(kt), whereas the space complexity remains O(k), the simulation has linear slowdown. 


With the use of the above lemma we can prove that every well formed QCA can efficiently be 
simulated by a well formed QTM, a question raised by Daniel Simon [51], John Watrous[57], and 
Christoph Dürr et al. [20, 21]. 


THEOREM 7.1 Every well formed QCA F can efficiently be simulated by a well formed QTM T. 
PROOF. For every QCA F there exists a QGCA G = ((0,1], n, M, S) which simulates F with linear 
slowdown (Lemma 6.10). Therefore, by Lemma 7.1, there exists a well formed QTM T' which 
simulates Fi with space-complexity O(k) and time-complexity O(kt). 


7.2 A Universal Quantum Cellular Automaton 


By now we have proven that a QTM can simulate a QCA and vice versa. We also know that there 
exists a universal QTM. This leads us to the conclusion that there exists a QCA that can simulate 
any other well-formed QcA. Still this is not the end of our investigations. The disadvantage of 
the above described construction is that we do not use the parallel computational power of a QCA 
when it has to simulate another parallel QCA. This is obviously not an efficient way of simulating. 
Here it will be shown that it is indeed possible to construct a QCA which can mimic the behavior 
of any other QCA in parallel. 

Our approach will be to use the known results about quantum circuits and universal quantum 
gates to construct a QGCA that we can ‘program’ with the input to simulate other possible QGCA. 
Because there are uncountable many QGCA we have to be satisfied with a simulation within a 
bounded error. This resembles the situation with the universal quantum gate. We will use this 
universal two bit gate U to approximate the M gate of a QGCA ({0,1},n,M,S). The resulting 
circuit after this transformation will be periodic both in the space and time direction. The universal 
QCA U will be a QCA which can be programmed to simulate any such circuit. 


7.2.1 Periodic Universal Gate Arrays 


Our aim is to simulate a QGCA G = (10, 1], n, M, S). For the simulation of the quantum gate M 
we will a universal two bit gate U. A circuit that only consists of U gates will be called a U gate 
array. The first transformational step will be to replace every M gate with a U gate array which 
approximates the initial gate M. Because of the symmetrical structure of QGCA, the resulting U 
gate array will be periodic in both space and time dimensions. Such a circuit will there be called 
a periodic U gate array. See Figure 7.1 for an example. In order to translate the ‘wiring’ of this 
array, another model will be introduced. 


7.2.2 Simulating the Wiring 


We proceed by defining two gates J and X that will mimic the wiring of a periodic U gate array. 
The gate I corresponds with the two-qubit identity gate whereas X describes a cross-over: 


Iz,y) = |x,y) and X|z,y) = ly c) 
for every x,y € {0,1}. By combining the 7, X and U gates in a wall-like, periodic structure, we 


can mimic any periodic U gate array. An of example such an periodic I XU array of this is shown 
in Figure 7.2. To describe these structures we use the following definition. 


DEFINITION 7.1 (periodic [XU array) A periodic IXU array F is defined by (m,n, A5, Bj) 
with A5, Bi € {I, X,U} for every j € Zm and i € Zp. 
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= U gate \> wire 


Figure 7.1: Part of a periodic U gate array which simulates a QGCA. Each box contains a U gate 
array that simulates a gate of the QGCA. The part shown here contains two layers of three original 
gates. The periodic pattern of the resulting U gate array is apparent. 


(| D I gate om X gate = U gate 


Figure 7.2: Example of a periodic XU array. The indicated areas correspond to the boxes of the 
periodic U gate array of Figure 7.1. This figure shows the simulation of one layer of three gates 
of the original QGCA. 
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A periodic [XU array is build by identical boxes of [XU arrays. The number m indicates the 
horizontal size (in terms of gates) of this box whereas n defines the number of double layers in 
each box. Each double layer consists of an A layer and a B layer. The function P5,,, describes 
the behavior of a row Z; of k boxes which acts on 2mk bits and and uses 2n time steps (measured 
in individual layers). This function is defined by: 


0 
Fomk = II st g Q Bi toam oe Q A$ modm 


i 1 jEZmk JEZmk 


EXAMPLE 7.1 The periodic I XU array of Figure 7.2 is defined by (2,2, A, B5) with 


Al =X ALU: 
B =U B? =] 
A= xX Ai =U 
Bead Bi=X 


Without going into full detail we can state the following lemma. 


LEMMA 7.2 For every QGCA F = ({0,1},n,M,S) there exists a periodic IXU array G = 
(n,m, A5, Bi) (with ged(m,n) = 1) which simulates F. This simulation is within an arbitrary 
small error. 

PRoor. Use the above described transformation from QGCA to periodic U gate arrays to periodic 
IXU arrays. The gcd(m,n) = 1 constraint is satisfied by adding a sufficient number of identity- 
layers which only consist of J-gates, thereby increasing the value of n until the desired number is 
reached. 


This lemma tells us that for every k € Nt the function Go, will simulate Fẹ. If € indicates the 
allowed error for simulating one M gate, the simulation of Fẹ will have an error of ke. Because 
€ can be made arbitrary small, ke can be made arbitrary small for a fixed value of k. The linear 
slowdown is expressed by the ‘vertical size’ of each block which equals 2n. If the simulation of the 
M gate requires some additional work bits for the U gate array, the space complexity will grow 
with a factor proportional to k. The reason for the gcd(m, n) = 1 constraint will become evident 
in the proof in Appendix D. 

If we express the cell values of the periodic [XU array at time t by |x, yb,- vu, po Uem a): 
the general behavior of the XU array is calibrated by the equations: 


A meam a 133) = le wj ) 7.1 
Bymodml¥j orja) = l ea) (7.2) 


for every i € N and j € Zkm. We will use this ‘calibration’ to prove that the universal automaton 
U embeds the same behavior. 


7.2.3 The Universal Quantum Cellular Automaton 


The universal quantum cellular automaton U is described by the following definition (an explana- 
tion of its semantics will be given below): 


DEFINITION 7.2 (U automaton) The U automaton is a QGCA (Q,n, M, P) with: 


state set Q = {0,1,<, DL, T; X VU" and n = 4 
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quantum gate M : Hoa — Has that obeys: 


M| 4, AY, ) RT , AS ) M| nn NT ) E | mn TA, ) 
M ? LA a) = |x, PAS ) M| , TAL ) = |<, FAN, ) 
Mx, TAN ) = |<, x, TAI, ) M T, TAT y) = |x’, AS, y) 
for every x,y € {0,1} and for every A € (1, X,U}: 
Alx,y) = |2",y’) 
neighborhood-scheme P = (n,o, ġ) defined by 
a(0) 23 c(1) 21 o(2) 22 a(3) 20 
$(0) = -1 o(1) =0 (2) =0 o(3) =1 


If we supply this automaton with the appropriate input, it will simulate a periodic [XU array. 
This is shown in the following analysis. 
For every N, K € N* the input: 


00 04,71 2 3 ip) 2N-1 
| wm, yO D, 00740, 4850, [HOO]... m 
. K 


with m € (TS X U^) and x,y € {0,1}, will evolve after two layers of the automaton Ufsw K 
to the configuration: 


0 1.1 ir3)2N-2 1 2N-1 
| TJ, dy jT; | TJ ees oe HT] 
JEL Kc 


with II9|z9, y9) = |x}, x£}) (the gate II is determined by 7 ="H97). At this point the z-values 
will travel step-by-step to the left, while the y-values have moved one place to the right and are 
now stationary (this is established by the “<”-symbol). 

The next important situation occurs after O(N) layers when the x and y values ‘meet again’ 
at the 7} gates, which are then to be simulated. This will result in the configuration: 


2r4..0 1 2 2 ip3]2N-1 
Lz Tj, OT 7O, Xv; Wy 3 TJ hes 
JEL 


with 11519525: = 1955954) and r} = T TEX 

By now the y-values go left in order to collide at the 73 gates with the z-variables. This process 
will be repeated in the obvious way. 

A careful examination of the behavior of U2) Kx on this input reveals that the overall time 
evolution is calibrated by the equations: 


2imod2N 2i 2i E 2i+1 ,2i4l 
MPPE [moa Kk 07i UJ) m Iz; VT) (7.3) 

2i+1mod2N 2i+1 |2i41 = +2 | 2i42 
IL, 2:2 | mod K Yo+i EJ it) = EAE TNTEY, (7.4) 


for every i € N, J € Zg, and v; = "IL". 
By using the calibrations 7.1, 7.2, 7.3, and 7.4, we can prove the following lemma. 


LEMMA 7.3 Every periodic IXU array F = (m,n, A5, Bj) with ged(m, n) = 1 can be simulated 
by the QGCA U. The time-complexity of the simulation is linear and does not depend on the size 
k of the periodic 1 XU array Foy. 

ProoF. If a ‘mapping’ between the four calibrations is possible, it follows that the two automata 
have the same evolution. Because gcd(m,n) = 1, there exists a positive integer a such that 
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an mod m = 1. With N = an and K = km every Fom,x can be simulated by Yan x. This is shown 
by the existence of a mapping which satisfies (assuming J = j — i): 


imodn 2imod2N imodn = [zit imod2n 


jmodm — j-i+ | 2$ | mod K jmodm j—-it+| 2233 | moa K 


for every i € N and j € Z. See Appendix D for the existence proof of such a mapping. Because a 
does not depend on k, the time-complexity is bounded by O(N?) = O(i) for the simulation of 2i 
layers of F. 


With this lemma we have reached the final conclusion of this thesis. 


THEOREM 7.2 There exists a proper universal quantum cellular automaton U which can simulate 
any well-formed QCA F with a bounded error. The extra costs of the simulation of Fj has time- 
complexity O(t) and space-complexity O(k). 

ProoF. This follows from the results of Lemma 6.10, Lemma 7.2 and Lemma 7.3. 


7.3 Conclusions 


We have defined a class of well-formed quantum cellular automata which are one-dimensional and 
circular bounded. It has been shown that every proper QCA resembles a periodic quantum gate 
array. This assures us that the global unitarity of a QCA can be reduced to the local unitarity 
of the gates of a QGCA. There exists a QCA U which can simulate any such array (by using a 
universal two-qubit gate U). The time-complexity of this simulation does not depend on the size 
k of the initial QCA. This proves that U is a universal quantum cellular automaton. 

The existence of a classical universal cellular automata was shown by Jürgen Albert et al. [1] 
and Bruno Martin [37]. Both constructions did applied to irreversible CA. Because the classical 
reversible cellular automata (RCA) are a subset of QCA, we know that every RCA can be simulated 
by the non-classical U automaton. If we want a universal RCA, we have to adapt Definition 
7.2 in such a way that it uses a universal three-bit gate (the Toffoli-gate for example). This 
is because there does not exist a classical two-bit gate U that is universal in its computational 
power[3, 4, 16, 18]. The author is unknown of an earlier construction of such an automaton. 

Theorem 7.2 shows that the model of QCA is a useful model of parallel quantum computing 
on its own. Because of their uniform structure, cellular automata are particular useful to bridge 
the gap between physics and computational theory. From this intermediate standpoint there are 
at least two directions possible, each with its own merits. By going from QCA to physics one 
obtains a powerful model to describe quantum mechanical systems. Recent work by David Meyer 
[39, 40, 41, 42] goes along this pathway. 

The other direction is to try to actually construct a controllable QCA. Several authors have 
suggested that QCA-like systems are more likely to be build than quantum Turing machines oriented 
structures [8, 34, 35, 36]. If such a construction would indeed be possible in the future, we would 
have equipped ourselves with a new remarkable tool. A tool whose computational power we are 
just beginning to unravel [11, 22, 56]. 


Appendix A 


Unitary Transformations 


Because of their central role in quantum computing, we will summarize the properties of unitary 
transformations [13, 25, 33, 47]. A linear transformation F is determined by its state set S and 
the corresponding transformational matrix Mp € C°*%. We will only look at countable state sets. 


DEFINITION A.1 A transformation F on S is unitary if and only if its corresponding matrix Mr 
obeys: MLMr = Mr M, = 1s, with lg the identity on the set S. 
A.1 Finite dimensional Transformations 


If Mp is finite dimensional (e.g. S is finite) MLMr = ls always implies Mr MT. = ls. We 
therefore can simplify the above definition. 


DEFINITION A.2 A finite dimensional transformation F on S is unitary if and only if its corre- 
sponding matrix Mp obeys: M} Mp = ls, with lg the identity on the set S. 

A.2 Infinite dimensional Transformations 

If we take the the infinite transformation F(i) =i+ 1 on S = N we obtain the following ‘inverse’ 


Ft: N+ SN, with Ft (i+1) = i for every i € N. We therefore have: FloF = I but not FoF! = T 
(F o F'(0) is not defined). In matrix notation this is illustrated by: 


0 1 0 .- 00 0 sss 100 
; 001- jo 02 Be 010 
Mp:Mp = 000- |1010- = 001 

but: 

00 0 010 000 
: 100 001 0 1 0 
Mr:Mp = 0 10 00 0 = 001 


A closer look at F reveals that F could not be unitary because F is not surjective. This is a 
necessary and sufficient restriction for infinite transformations with MLM. F = ls to be unitary: 


DEFINITION A.3 An infinite dimensional transformation F on S is unitary if and only if its 


corresponding matrix Mp obeys: MLMr = ls (with 1g the identity on the set S) and F is 
surjective. 
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A.3 Exponential Expressions 


DEFINITION A.4 Let H be a square matrix, by e" we mean I + H + iH? +.--. Therefore for 
every t € C: 
1 
= — 5 k . k 
exp(t. H) = 22 t°- H 
k=0 


EXAMPLE A.1 
a a a oa 


DEFINITION A.5 A matrix A is called hermitian if and only if A! = A. 
Finite unitary and hermitian matrix correspond to each other in the following way. 
1. For every hermitian matrix H, exp(iH) is a unitary matrix. 


2. For every unitary matrix U, there exists a hermitian matrix H such that: U = exp(iH) (in 
general this H is not unique). 


For this reason the analogy of hermitian matrices to real numbers and unitary matrices to the 
complex numbers with norm one is made (except for commutativity of multiplication). 


Appendix B 


Proving Well-Formedness 


In this appendix we will prove that if the general function Fz of a normalized QCA is not well- 
formed, there exists a K € N+ such that Fẹ is not unitary. This result is used in Lemma 4.3. 


LEMMA B.1 Let F = (Q, N, f) a normalized QCA. If there exist two values X 4 Y € Q? such 
that 


(Fz(X),Fz(Y. # 0 


for every i € Z, then there exists ak € Nt with Fẹ not unitary. 
PROOF. We use zx; and y; € Ha defined by 


Fz(X) = Fz (® x) - 9n; and Fz(Y) = Fz le x) =u 
ieZ i€Z i€Z i€Z 
with X;,Y; € Q and Xo Æ Yo such that 
for every i €Z: (xi yi) #0 (B.1) 
We will now prove that, given these X and Y, we are able to construct ak € N+ and a X^, Y' € QF, 
with X' L Y' and Fy| X) £ Fi|Y^. 


Define r = |N| the size of the normalized neighborhood N = (0,1,...,r —1). Take a,b,c and 
d (with a <b <0 < c< d), such that: 


Xa:(a+r-1) = Xb:(b4r—1) Yas(at+r—1) ES Yb:(b+r—1) 


Xe(ctr—1) = Xd:(d4+r—1) Yo:(e+r—1) = Yas(d-+-r—1) 
Because there are ‘only’ |Q|?" different combinations of [Xi:(¢t+r—1)3 Yi:G+r—1)], this is always pos- 
sible for a > —|Q|?" and d € |Q|?". With these values we define: 

Xr = Xa:(b-1) XR = Xa(a-1) 

Yr = Yu Yr = Yo:(a-1) 


This can be visualized by: 


XL XR 
i 
X=.. Xa: 8D Xoi @QX, 0 X09 @---® Xe- Q Xq4 QX 
Y =... a8: 9Y,4 II- IYD -Q Y- Q-Q Yı QYq 
SS ——— 


Now we distinguish the following three possibilities: 
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1. Xr # Yr: Take X' = Xz, Y’ = Yr both elements of Q^ with k = b — a (thus X’ L Y"). For 
0<i<r-—1I1 we have: 
Kage X; = Xati = Xoti 


and therefore: 


F(X’) 


(2s) 


jEZk iEZr 


- Qi (Q kam) 


jEZk iEZr 
= (Fz(X))a:(6-1) 


and for the same reasons also: F;(Y^) = (Fz(Y)),,4, ,). This leads us to: 


b—1 


(Fe(X"), F (Y) = (O00). (FeO axe) = IE Gov? 


i—a 


By (B.1) we now know that F(X’) L F(Y’), which shows that Fẹ is not a unitary trans- 
formation (with k < |Q["^). 


2. Xr z Yr: With X’ = Xr, Y' = Yr and k = d — c, the same reasoning as with Xr, Æ YT, 
holds (again k < |Q["^). 


3. Xr, = Yr and Xr = Yr: First we will prove Xa (b+r-1) = Ya:(b+r—1)- We already know: 
Xa(p—1) = Yax(p-1). With induction (given 0 < i € r — 1 and X4(5,14 1) = Ya:(b+i-1)) 
it follows that Xy; Xari Yati Yi; holds. Likewise we can prove: X¢.(d4r—1) = 
Ya(apr-i) Now take X' = Xy(o4,1,; and Y' = Yg(e45.1) € Q^, with k = c4 r — b. 
Because b € 0 € c-r — 1 and Xo Z Yo, we have X’ L Y’. Below we will use: (because 
Xb:(b4r—1) = Yp.(b4r—1) and Xo:(c4+r—1) — Y cues): 


for every 0<i<r—1: Xl ppi = hee = Yopi = Yo ppi 


and also for every r < i < 2r — 2: 


/ / £ p E t EM f =. f 
X. p = Xkti-r T X; = Xpti-r — Yosi-r E Yi = Yi nios = Y. pn 


To summarize: X’ muy ;if0< i< 2r —2. This gives us: 


c—b4i c—b-4ri 
EU) = @1(@x i 
j€Zx icZ 


c—b—1 k—1 
= f (8x) 8 [ f (& xs 
j=0 i€Z. j=c—b icZ 


g (g=) (8,87) 


icZ 


If we combine this with: 
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the inner-product of F(X’) and F(Y’) equals: 


(FAX), FAY) = D (Bx) | D: (Dr) ) 
jab 


j=b i€Z, i€Zy. 


= ((F (X)u qe» Cz) 5) 


c—1 


= (zi, Yi) 

i=b 

Again, (B.1) shows us that F(X’) £ F(Y’), which shows that Fẹ is not a unitary trans- 
formation (with k < 2|Q["" +r). 


All possibilities are covered by 1, 2 and 3. This proves that there exists a k < 2|Q[^" +r, for which 
Fy, is not unitary: the QCA F is not well-formed as defined in (4.7). 


Appendix C 


Simulating Neighborhood Schemes 


In this appendix we give the proof of Lemma 6.6 about the simulation of arbitrary neighborhood 
schemes by periodic QGCA which only use the Shift neighborhood scheme. 


LEMMA C.1 For every QGCA F = (Q,n, M, P) there exists a periodic QGCA G defined by 
G = (Q,n, u, Mo, ..., Mu, $,..., 9) 


which is shift equivalent with F. 
PROOF. Given the neighborhood scheme P = (n,a,¢), we will use the labels: 


e B= (ó(i) |i € Zn} 

e r = max(®) — min(®) 

e Wi = (c(t) | t € Zn and ó(t) = min(9) +i} forevery0 iz r 
e m(P) = Y igi: vil 


(Therefore if m(P) = 0 we have yo = Zn and m(P) = 1 implies |yo| = n — 1 and |vi| = 1.) 
For every neighborhood scheme P = (n,c,$) and value c € Zn, we define the following standard 
procedure: 

Construct an alternative neighborhood scheme R = (n, s, v) by 


«G)-—c(i) and (i) = H(t) — 1 (2) 
Next, we define two ’permutation gates’ II: Zn > Zn and 7 : Zn > Zn by: 


n—1 ifi=c 


n) = c ifi2n-1 
i otherwise 
and (i) = II(; — 1) (such that z(II(;) + 1) = i for every i € Zn). We now have for every k € N*: 
II o S, ox o R, = (n, II(i), 0) o (nii -- 1, 15 1(2)) o (n,m(1),0) o (n, sli), p(i)) 
= (n,1(1),0) o (n,i - 1,1, 1(2)) o (n, co m(i),q om(1)) 
= (n,1(1),0) o (n,com(i +1), 15 1(2) + g om(i 4 1)) 
= (msom(II(i) + 1),1, 1 ° H(i) + p o TCi) + 1) 
= (n,o(i), (i) 
= P, 


This shows that the u + 1-periodic QGCA G’ = (Q,n,w+1,7,Mo oll,... , My, R,S,... , 8) is 
identical with the -periodic QGCA G = (Q,n, u, Mo,... ,My-1, P, $,..., S), and thus G = G’. 


65 


66 


APPENDIX C. SIMULATING NEIGHBORHOOD SCHEMES 


We will use this procedure to ‘transform’ the QGCA F into the desired periodic QGCA. For every 
-periodic QGCA G = (Q,n, u, Mo,...,My-1, P, $,... , 9), we distinguish the following (mutual 
exclusive) possibilities: 


1. 


{0} = vo and n = 1: The neighborhood scheme P = (n,o,¢) has o(0) = 0 and ¢(0) = t 
(with t € Z). We therefore have S o S*~! = Pj for every k € N+. This proves the j-periodic 
acca G' = (Q,n, u, Mo,... , My 1, $,... , S) to be shift-equivalent with the initial QGCA 
G; in other words G' = G. 


. {0} = vo and n > 1: Take c such that o(c) = max(W) and apply the procedure. The new 


neighborhood scheme R has m(R) = m(P) — 1. 


. {0} € ve: Take c with o(c) # 0 and ¢(c) = min(W) and apply the procedure to get 


the neighborhood scheme R. Now we have yè = (c(c)) = {a(c)} Z {0} combined with 
m(R) =m(P)+n-1. 


. {0} Z vo and m(P) > 1: Take c such that o(c) = max(y,) and apply the procedure. For 


the neighborhood scheme R we have m(R) = m(P) — 1 and (again) (0) Z yè. 


. {0} Z yo and m(P) = 1: This implies: Yo = (1,...,n — 1}; v1 = {0}; 9 = (tt 1}; 


olc) = t+ 1; o(c) = 0, with t € Z and c € Zn (therefore (i) = t + 1e(i)). Define the 
permutation-gate II(7) = c (i) — 1. It now holds for every k € Nt: 


Io Spo St” = (n,I(i),0)o (n,i-- 1,14 1(i)) o (n, i,t) 
= (nII(2) + 1, 15 4 (II(2))) o (n, i,t) 
= (n,o(t), 1-(4) +t) 
= B, 


This shows that if we define a u-periodic QGCA 
G' = (Q,n, p, Mo 2 TL M3, Mua 9,- 5) 


we have finally reached the shift-equivalence G = G”. 


If we start with the initial QGCA F and apply the above transformation recursively, we will always 
reach (after a finite amount of steps) possibility 5. This means that we have constructed the 


desired -periodic QGCA G, with F = G. 


Appendix D 


Mapping the Calibrations 


We will prove the existence of a mapping as used in Lemma 7.3. 


LEMMA D.1 For every m,n € N* with gcd(m, n) = 1, there exists an a € N* and a mapping II: 
of A; and B: which respects the equations (with N = an and K = km) 


imodn = I[24mod2N imodn __ [zit imod2nN 


jmodm j-i+ | 25 |modK jmodm — j—i+ (24! | mod K 


for every i,k € N and j € Z. 
PROOF. In order to ensure a correct mapping between the A values and II, we have to prove that 


if: 
2imod2N = 2pmod2N 
TIE mod K = q p| mod K 
then: 
imodn = pmodn 
jmodm = qmodm 


for every i, p € N and j,q € Z. 
By K = km and defining a such that N = an = 1 mod m, the initial conditions can be 
restated to: 


2(;— p) mod 2an = 0 
; F 20 2p 
(j—4)—- (i— p) + ES - E mod km = 0 
2an 2an 


The first condition leads to p = i+ Aan (with A € Z). As a result the second condition becomes: 
(j—q)-- A(an — 1) mod km. = 0 


Therefore the the restrictions (i — p) mod n = 0 and (j — q) mod m = 0 are satisfied for every 
i,p € N and j,q € Z. 

'The same reasoning holds if we want to prove the existence of a consistent mapping between 
the B values and II. 
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Appendix E 


Sources of Information 


The writing of this thesis has benefited greatly from the information which is nowadays available on 
the internet. The pace of the recent developments in quantum computing is unthinkable without 
the existence of archives such as quant-ph or qcomp. What follows is a personal selection of the 
sites and home-pages which I found particular useful. 


http: //xxx.lanl.gov/archives/quant-ph 
The e-print archive of the Los Alamos National Laboratory has shown to be the center of 
publications regarding quantum computation. 


http://babbage.sissa.it/quant-ph/ 
'The European mirror-site of quant-ph. 


http://feynman.stanford.edu/qcomp/ 
The Quantum Computation Archive at the Stanford University provides a large but well- 
selected database of articles solely about quantum computation. 


http://eve.physics.ox.ac.uk/QChome.html 
'The home-page of the Quantum Information Group at the University of Oxford. With on- 
line elementary introductions about quantum cryptography, cryptoanalysis, computation 
and communication by Artur Ekert and David Deutsch. 


http://vesta.physics.ucla.edu/^smolin/index.html 
'The colorful Quantum Information Page of John Smolin at University of California, Los An- 
geles. 


http://www.cwi.nl/^berthiau/ 
André Berthiaume's home-page with his introduction in quantum computation [7]. 


http://www.iro.umontreal.ca/labs/theorique/index en.html 
The Laboratory for Theoretical & Quantum Computing at the Université de Montréal. Gilles 
Brassard's extensive Bibliography of Quantum Cryptography can be found at the English 
home-page of Claude Crépeau. 


http://www.research.ibm.com/quantuminfo/ 
The research group at IBM with an on-line article about quantum teleportation. 


http://www.physik.uni-ulm.de/^sam/home.html 
Starting point of the tutorial on quantum computation by Samuel Braunstein. 


http://publish.aps.org/PRA/prahome.html 
Site of Physical Review A, the journal where the field of quantum computation has found 
its niche. 
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http://alife.santafe.edu/alife/topics/cas/ca-faq/ca-faq.html 
Frequently Asked Questions about cellular automata; edited by Howard Gutowitz. Keeps 
also track of the ongoing discussions in the newsgroup comp.theory.cell-automata. 


http://cs.ua.edu/graduate/lusth/qca/ 
The Quantum Coupled Architecture Home Page which deals with a number of subjects 
closely related to QCA. 


http://aerodec.anu.edu.au/~qc/index. html 
The Quantum Computing home page of the Australian National University, including a list 
of paper reviews. 


http://altavista.digital.com/ 
“Tf all else fails, try consulting Alta Vista.” Current statistics (Summer 1996): 
e “cellular automata”: 5037 documents matching the query 
e “quantum computation”: 990 documents matching the query 


e “quantum cellular automata”: 40 documents matching the query 
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